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UNSTEADY FLOW CONDUITS WITH SIMPLE SURGE TANKS 


ABSTRACT 


The basic principles the operation simple surge tanks have been 
known for decades. American and European literature contains many articles 
devoted this subject. the fact that the fundamental flow equa- 
tions are simple concept, they are non-linear and hence difficult solu- 
tion. The bulk the existing literature with one relatively recent exception(1) 
has been concerned with numerical and graphical methods solution. This 
paper presents general solution the simple tank problem terms 
single dimensionless parameter for the case instantaneous and complete 
load rejection. The mechanical differential analyzer the University 
California was used obtaining the solution. Seventeen curves covering the 
range from the case practically damping that almost critical damp- 
ing are presented. Whereas nothing new principle presented, be- 
lieved that the remarks below constitute refinement existing solutions 
and present numerical data covering wide range conditions concise 
manner. 


INTRODUCTION 


The problem under consideration involves the prediction the water sur- 
face elevation the surge tank, depicted schematically Fig. func- 
tion time following rapid and complete cessation flow the conduit. 


Note: Discussion open until July 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1933 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. February, 1959. 


Prof. Eng., Univ. California, Los Angeles, Calif. 

Asst. Prof. Eng., Los Angeles State College, formerly Mechanical Engr., 
Los Angeles Dept. Water and Power, Los Angeles, Calif. 

Computer staff Beckman/Berkeley Division, formerly Graduate Research 
Asst., Univ. California, Los Angeles, Calif. 
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Theory 


The fundamental equation unsteady flow the conduit, based the so- 
called “rigid water column” analysis, is, consideration Fig. 


(1) 


which measured positively upward. Very simply interpreted, the equa- 
tion states that the driving head less the friction head equal the inertia 
the water column the conduit. assumed that the inertia and friction 
the water column the tank are negligible comparison with those the 
pipe. 

The continuity condition the conduit-tank junction requires that 


(3) 


for the case instantaneous and complete closure. 
Differentiating (3): 


(4) 


and are eliminated between (3), (4) and (2) upon arrangement terms: 


simplifying: 


account for damping during flow reversal the conduit necessary 
write eq. (5) the following form. 


2 


This merely precludes the possibility the damping term having opposite 
sign from the other terms which case energy would appear added 
the system dissipative forces. 

Equation (6) may rendered dimensionless the following changes 
variables: 


i 
: 
{ 


Let 


worth noting that the period undamped oscillations the conduit- 
tank system. 
Further let 


(b) 


Equation (6) now becomes 


Multiplying 


The dimensionless coefficient further designated the symbol 
the damping factor for the normalized equation. Finally 


The initial conditions, corresponding the instant the valve closed are: 


the 
which the driving head for steady flow. 


Computer Solutions 


diagram the differential analyzer mechanization shown Fig. 
Our remarks will restricted this particular problem since the theory 
the differential analyzer has been presented the literature.(2,3) For the 
most part, standard techniques were used obtaining solutions Eq. (8). 
The equation was written the form 


and the right-hand side this equation constructed repeated integration, 
constant multiplication, and summation. The use integrator (No. 
multiply constant standard where large number parameter values 
are used. clutch “disconnect” between integrators and order 


was the only special feature involved. Seventeen solutions eq. (8) 


became negative (so obtain 


dimensionless plot indicating the height and time the first maximum surge 


CONCLUDING REMARKS 


was found possible express the differential equation describing the un- 
steady flow surge tank system sudden outflow cessation terms 
dimensionless variables and single dimensionless “damping” parameter. 
This equation was solved mechanically and the results are presented 
form which makes possible the determination the maximum surge, the 
time its occurence and approximation the time complete cessation 
water pendulation, for most simple surge tank system configurations. 

The computer solutions were substantiated within accuracy numeri- 
cal analysis typical system.(4) experimental model study the 
University California, Los Angeles (unpublished) likewise indicated sub- 
stantiation the computer solutions the extent that the height and time 
the first maximum surge were found excellent agreement with the 
above theory. However, the rate damping observed actual flow has 
been tentatively found variance with that predicted the theory. 
addition, the conditions stability simple tank system were found ex- 
perimentally also disagreement with the Thoma criterion which 
widely quoted the literature.(5,6) These latter matters are currently sub- 
ject investigation. 
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SYNOPSIS 


There definite evidence that the precipitation falling over the surface 
the Great Lakes considerably less than the precipitation falling adjacent 
land areas. Evaporation computed the water budget method reduced 
the same amount the precipitation. The annual amount net evaporation 
from Lake Ontario computed the water budget and mass transfer methods 
the order 2.0 feet. 


INTRODUCTION 


The advent the St. Lawrence Seaway projects has stimulated the diverse 
interests landowners, shippers, economists, and engineers, all whom 
have been concerned about the levels the Great Lakes for some time. Yet, 
despite their attentive concern, many are not aware that the outflow the St. 
Lawrence River and, thus, the water levels Lake Ontario, have been com- 
pletely regulated the Iroquois Dam control structure since 1958. 

the responsibility the Corps Engineers, Army, ascertain 
for the United States that the regulation Lake Ontario provides optimum 
benefits for navigation, power, and riparian interests. this end, the 
Lake Survey has been studying techniques forecasting supplies the Great 
Lakes for future use effecting the most beneficial regulation plans. These 
studies determined that, all the hydrologic factors affecting the levels 
the Great Lakes, evaporation was the factor about which the least was known. 

Evaporation the process which water converted vapor. The 
molecules that make water are continuous motion. Some molecules 
water have sufficient momentum break through the surface and enter the 


Note: Discussion open until July 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1934 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. February, 1959. 

Presented the February 1958, ASCE Convention Chicago, 

Maj., Lake Survey, Corps Engrs., Detroit, Mich. 
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air. The rate which these water particles leave the water and enter the 
adjacent air depends upon the heat supply the water and the condition the 
air. Concurrently, molecules water vapor the air strike and enter the 
water; that is, they condense. hydrology, evaporation considered 
the net loss water mass. 

There are several methods computing evaporation. For the Great 
Lakes, the water budget method has been utilized most frequently. More 
recently, evaporation has been correlated with the vertical distribution 
moisture the air and the intensity turbulent mixing. This the so- 
called mass transfer theory. Another, called the energy budget method, 
takes into consideration solar radiation, atmospheric radiation, advected 
energy, sensible heat transfer, and other factors which affect the total 
amount heat available convert water into vapor. Computing evaporation 
the water budget and mass transfer methods will discussed herein. 


Water Budget Method 


The water budget method determining evaporation Lake Ontario 
simply matter solving the general hydrographic equation for evaporation: 


where 


the evaporation from the surface Lake Ontario during the inter- 
val 


the change the mean level Lake Ontario during the interval 


the runoff water from the land surface draining into Lake 
Ontario during the interval 


the precipitation falling directly the surface Lake Ontario 
during the interval 


the inflow Lake Ontario from Lake Erie, the flow the 
Niagara River plus the Welland Canal, the interval 


the outflow from Lake Ontario, the flow the St. Lawrence 


the sum all factors affecting the rise fall the level Lake 
Ontario not contained the other factors the equation. 


Although this method determining evaporation simple, the determina- 
tion the component factors which make the equation quite difficult. 
studying Lake Ontario, the amount hydrological and meteorological data 
available the period prior 1934 was limited, twenty-year period 
record from 1934 through 1953 was used. Evaporation was computed 
monthly units because the majority data was available monthly basis 
and because was intended that the results the study would used 
regulation plans for Lake Ontario where thirty-day forecast the water 
levels desired. 
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Normal changes levels Lake Ontario 


The Lake Survey, coordination with the Water Resources Branch, 
Department Northern Affairs and National Resources, Canada, has deter- 
mined Lake Ontario end-of-month levels. From these data, the average 
monthly changes Lake Ontario levels for the period 1934 1953 were 
computed and are tabulated Column Table 

the period record, 1934-1953, the all-time mean annual low (1860- 
1957) 243.54 feet was established 1935, and the all-time annual high 
(1860-1957) 248.01 feet was established 1952. During this maximum 
secular climb Lake Ontario, the average lake levels rose during the period 
January through May and fell during the period July through December; the 


TABLE 


Evaporation Feet Computed Water Budget Method 
Period Record, 1934-1953 


JANUARY 


1-2 


MONTH 


FEBRUARY 


MARCH 


APRIL 


MAY 


JUNE 


JULY 


AUGUST 


SEPTEMBER 


OCTOBER 


NOVEMBER 


DECEMBER 


| 


maximum monthly rises occurred March and April, and the maximum 
monthly declines occurred August and September. 


Runoff 


The runoff precipitation from the drainage basin Lake Ontario 
extremely important factor analyzing lake level fluctuation. Preliminary 
studies have shown that there small lag from the time precipitation 
until the runoff reaches the lake. Although several the tributary basins 
have great number small lakes and the basin the Trent River par- 
tially regulated, there does not appear any marked degree storage 
the Lake Ontario basin. the 34,795 square miles the basin, the lake sur- 
face consists 7,520 square miles and the land drainage area consists 
27,275 square miles; therefore, not difficult visualize why runoff 
important. 

The majority runoff from the land area into Lake Ontario the form 
stream flow from the major rivers tributary the lake. The files 
United States and Canadian agencies were examined for available stream flow 
records. Usable records for nine major rivers which include 16,500 square 

miles per cent the Lake Ontario land area (excluding that area tribu- 
tary the Niagara River above Niagara Falls) were found. Drainage areas 
with available records were incorporated into the near drainage areas with 
surface records, see Figure 
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Figure 1. Drainage Areas, Lake Ontario 


The areas for which gaged records were not available were divided almost 
evenly between upland and coastal regions. Although there not much 


fir 
TRENT RIVER 
CREDIT RIVER 
LAKE ONTARIO 


variation the topography the Lake Ontario Basin, except the eastern 
portion where mountain foothills are located, care was taken incorporate 
ungaged areas with the gaged areas that had the same general topographical 
characteristics. For instance, the ungaged area the coastal region 
northern Lake Ontario was incorporated with Duffin Creek which also lies 
within the coastal region. The ungaged upland regions east the Moira 
River were incorporated with the Moira River basin. The size the gaged 
drainage areas the nine rivers and the areas incorporated with the records 
these rivers are given Table The Credit and Humber Rivers and Duf- 

fin Creek had limited records and the statistics flow were mathematically 

extended the period 1934-1953. 

The average monthly runoff for the lake basin measured feet Lake 
Ontario for the period 1934-1952 given Column Table 
Although there are certain some errors measuring the total run- 
off, the gaged runoff the most accurately measured variable the hydrolog- 
cycle for, instead being point measurement within large area (such 
the measurement precipitation), the gaged runoff effectively integrates 

the entire area from which the measured flow originates. 

Runoff complex phenomena which depends upon many variables. The 
observed monthly runoff for the Lake Ontario Basin function the supply 
water available for runoff and the losses which prevent the supply from 
reaching the lake. The supply water available during month directly 
related the precipitation during that month, the interflow, ground water, 
and channel storage, and the water equivalent the snowpack; any. The 
losses are due interception, depression storage, evapotranspiration, and 
replenishment soil moisture. 

Indexes these most important variables affecting runoff were deter- 
mined and mathematical correlation historical record with these indexes 
were obtained which gave results rather high correlation.(1) 

The relation runoff precipitation proved useful tool gaining 
insight into the runoff from the basin. This relationship for each the 
drainage areas and for the entire basin given Table Some very inter- 
esting results are shown. One, the extremely large runoff factors for the 
months March and April. This demonstrates conclusively that knowledge 
the storage precipitation the form snow and the subsequent snow 
melt vital importance complete understanding spring runoff 
obtained. Another, the magnitude the yearly runoff factor for the 
Lake Ontario basin. The value 48.0 per cent appreciably higher than the 


value 30.7 per cent obtained John Freeman his study, “Regulation 
Great Lakes.” 


Precipitation over Lake.Ontario 


The amount precipitation over the lake affects the lake level immediately. 
Heretofore, the weighted precipitation catches recorded selected weather 
stations the perimeter the lake have been regarded the precipitation 
falling over the lake. The perimeter precipitation shows marked reduction 
from the precipitation falling over the land drainage area, yet seems that, 
due the cooling effect the lake the summer and the fact that there are 
diurnal heating orographic effects, the precipitation falling over the lake 
much less than indicated the perimeter stations. 
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late years there have been several concerted attempts determine 
the depth precipitation falling over large bodies water the same 
order magnitude the depth falling over the adjacent surrounding land 
areas. the Great Lakes area, precipitation has been recorded St. James, 
Beaver Island, northeastern Lake Michigan, for twenty-five years during 
the period 1911-1956. comparison the precipitation falling Beaver 
Island with the precipitation falling the adjacent land areas indicates that 
not only the depth precipitation less over Beaver Island, but that there 
seasonal variation. 


Northeastern Lake Michigan studies 


investigate the theory that the depth precipitation falling over large 
bodies water less than the depth falling over adjacent surrounding 
land areas, study precipitation over northeastern Lake Michigan was 
instituted the United States Lake Survey. The data for the study were ob- 
tained jointly with the Weather Bureau and with the cooperation the 
Coast Guard. Storage precipitation gages were installed six islands 
during November 1952, the Weather Bureau, see Figure These 
special island storage gages were inspected and the catch was measured 
semiannually. The results the catch the special island gages have been 
compared with network ten shore control stations. 

The topography the six islands varies tremendously and the ground 
cover and exposure the gage sites also differ greatly; yet, the catch 


LAKE SUPERIOR 


BEAVER 


SS 


MamTOU 
8. MANITOU 


Legend: 


PRECIPITATION 
STORAGE GAGES 


Figure 2. Precipitation Gage Locations, Northeastern Lake Michigan 
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five the six island stations for the period November 1952 November 
1957, compares very favorably with that St. James, Beaver Island. The 
average precipitation the five large island stations for the five year period 
103 per cent the precipitation recorded St. James. However, the 
catch Ile aux Galets, the smallest the six special island stations with 

area 2.9 acres, was appreciably less than the catch the other island 
stations. aux Galets very small flat sand spit and completely de- 
void vegetation, whereas the other island stations are fairly heavily vege- 
tated and relatively large area. 

The greatest difference precipitation between Ile aux Galets and the 
other island stations the winter period. This believed caused 
the extreme exposure the gage, which reduces the amount catch 
snow. 

The records during the summer period show that the depth precipitation 
Ile aux Galets approximately the same for South Fox Island. The 
average precipitation Ile aux Galets and South Fox Island during the period 
May through October, 1952-1957, when comparable data are available, 
per cent that recorded St. James. This reasonable, because Beaver 
Island large enough for definite diurnal heating occur and the convective 
currents generated could cause additional precipitation compared with the 
smaller islands. 

Nevertheless, determine the size and vegetative cover the island 
had definite effect the amount precipitation, additional gage was 
installed October, 1956, Shoe Island, another sand spit 2.3 acres size. 
Figure shows the gage location Shoe Island. The results the first year 
operation corroborate the reduced catch aux Galets, indicating that 
the rainfall over the small islands, and consequently over the lake, re- 
duced from that recorded Beaver Island. 

fully realized that five years and less are very brief periods which 
base specific conclusions regarding precipitation. But, general trends can 
ascertained and approximate relationships can determined. the 
other hand, the twenty-five years record St. James are much more reli- 
able. The ratio precipitation recorded St. James that recorded over 
the surrounding land area given Column Table These ratios show 
definite reduction rainfall, with the greatest reduction occurring the 
summer months. This what may expected, because there are fewer 
frontal storms the summer months and those that pass over the lake are 
dampened due the cooling effect the lake surfaces. The convective type 
precipitation pronounced over land during the summer period. the win- 
ter, there are more frequent passages frontal storms, dampening effect 
the lake, and squall-type storms. Therefore, one would expect the 
precipitation over the lake approximate that occurring over the land. The 
difference the air and water temperature the vicinity Beaver Island 
illustrated Figure can seen that the average monthly air tem- 
perature higher than the average monthly water temperature during the 
period April through August, accounting part for the reduced rainfall over 
the lake during this period. 

previous paragraph was stated that the precipitation over the lake 
was indicated less than that recorded St. James. From the records 
South Fox Island, [le aux Galets, and Shoe Island, seemed reasonable and 
conservative assume the reductions shown Column Table Thus, 
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Figure Precipitation Gage, Shoe Island 


was assumed that the product the two ratios, the ratio St. James land 
area and the ratio small island stations St. James, would give approxi- 
mate solution the reduction precipitation falling over northeastern Lake 
Michigan. plot smoothed curve drawn through the values the product 
the two ratios, the ratio precipitation over water precipitation over 
land, appears Figure and the “smoothed” monthly values used reduce 
land precipitation over northeastern Lake Michigan are given Column 
Table 

his paper entitled, “The Influence Lakes and Urban Areas Radar 
Observed Precipitation Echoes,” Bulletin American Meteorological 
Society, Pearson gives the conclusion that, “Observation radar film 
echoes the vicinity over Lake Michigan revealed evidence that the 
lake, general, discourages the formation air mass showers.” This fact 
definitely evidenced the outline radar echoes depicted Figure 
addition, data collected the Jacksonville District, Corps Engineers, 
precipitation stations Lake Okeechobee, the second largest fresh lake 
the United States with area 730 square miles, reveal reduction 
rainfall over the lake compared with the catches the perimeter the lake. 
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TABLE 


Precipitation Ratios, Northeastern Lake Michigan 


Islands 
Surrounding Land 


St. James 
Surrounding 
Land 


Small Islands 
Surrounding 
Land 


Small Islands 
St. James 


January 
February 
March 
April 


fax 


September 
October 


November 


December 


Yearly 
Average 


*Hypothetical 


The Weather Bureau office Key West has approximated their rainfall 

per cent that falling the Florida mainland. Thus, there conclu- 

evidence that the precipitation falling over large bodies water re- 

One may reasonably ask the ratios obtained for northeastern Lake 
Michigan are applicable Lake Ontario. Evidence indicated that approxi- 
mately the same ratios are applicable. Generally, the same prevailing winds 
hold for the whole Great Lakes region. comparison precipitation the 
perimeter stations with the land drainage area stations Lake Michigan and 


4 
Month 
| 
1 
| q 
4 


MONTHS 


Legend: 


AVG. LAKE ONTARIO 
AVG. LAKE MICHIGAN, VICINITY JAMES 


Figure Average Monthly Air-Water Temperature Differentials 
Degrees Fahrenheit, Lakes Michigan and Ontario 


Lake Ontario showed that the ratios are remarkably the same. shown 
Figure the reduction precipitation over both lakes follows the same pat- 
tern. Furthermore, comparison the average monthly water temperature 
minus the average monthly air temperature revealed that almost exactly the 
same pattern followed both Lakes Michigan and Ontario, see Figure 
above. This indicated that the lakes discourage the formation air mass 
showers during the same period the year. Accordingly, factors developed 
for reducing the amount precipitation falling over northeastern Lake 
Michigan will apply Lake Ontario. 

The average monthly precipitation falling over the Lake Ontario land 
drainage area constant, varying but slightly from the average monthly 
mean 0.24 foot. The average monthly over-land precipitation shown 
Column Table above. This average monthly precipitation over land, 
multiplied the ratios previously computed for reducing the precipitation 
over land the precipitation over water Column Table gave the aver- 


age monthly depth precipitation falling Lake Ontario, see Column 
Table 


Inflow Lake Ontario from Lake Erie 


The inflow Lake Ontario from Lake Erie the flow the Niagara 
River plus the diversion water through the Welland Canal Canada. The 
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Figure Ratios, Precipitation Over Water Precipitation Over Land 


Welland Canal diverts water for navigation and power from Lake Erie Port 
Colborne and returns Lake Ontario Port Weller and Port Dalhousie. 
The Niagara River flows used this study were computed the new Morrison 
Street rating curve. The Morrison Street gage located the Maid the 
Mist Pool below Niagara Falls and was the only gage the Niagara River 
which was comparatively free from ice and weed effects during the period 
record. interesting that the flows used this study differ from the pub- 
lished flows the Niagara River which occurred high Lake Erie levels 
much 10,000 cfs per month. Revision the published Niagara River 
flows being studied the United States Lake Survey. 


Outflow the St. Lawrence River 


The Coordinating Committee Great Lakes Basic Hydraulic and 
Hydrologic Data(2) has accepted the premise that the Rapide Plat weir the 
best available meter the flow the St. Lawrence River and that the flow 
over the weir can correctly measured utilizing the stage the Lock 
gage. The committee, therefore, adopted the Lock stage-discharge rela- 
tionship the basic relationship for the determination the river flow. The 
Lake Ontario outflows utilized this study are those adopted the Coordi- 
nating Committee Great Lakes Basic Hydraulic and Hydrologic Data. 
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Figure Outlines Precipitation Echoes Taken 30-Minute 
Intervals. Courtesy Illinois State Water Survey 
Division, Urbana, Illinois 


The average monthly inflows Lake Ontario from Lake Erie minus the 
average monthly outflows the St. Lawrence River, feet Lake Ontario 
levels, for the period 1934-1953, are shown Column Table above. 
These values not vary greatly and there excellent correlation between 
the inflows and outflows when they are analyzed statistically. 


Other factors 


All other factors which affect the amount evaporation computed 
water budget equation not previously discussed were included the item 
The most important these factors were the ground water flow and 
the thermal expansion water. 
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Figure Ratio, Precipitation Lake Perimeter Stations 
Precipitation Over Land, Lakes Michigan and Ontario 


Ground water flow 


The average level Lake Erie 572.3 feet, whereas, the average level 
Lake Ontario 246.0 feet. was possible that there might considerable 
underground flow between the two lakes due the very large potential head. 
This matter was discussed with members the Geological Survey and 
the consensus was that there was not appreciable ground water flow into 
Lake Ontario. 


Thermal expansion water 


mass water heated and cooled, will expand and contract. The 
greater the depth, the larger the thermal expansion and contraction. The 
average depth Lake Ontario 264 feet. However, already shown, 
definite thermocline exists Lake Ontario. This thermocline varies depth 
both with time and position, but the maximum depth does not exceed feet, 
below which the water generally temperature 39° The maximum 
monthly average water temperature difference between two consecutive 
months was 11° see Table above. The lake expands vertical direc- 
tion only, and the expansion function the change average density 
the water above the thermocline. The maximum amount expansion and con- 
traction about .01 foot per month. Consequently, the thermal expansion 
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water was not considered the water budget equation. 

addition the ground water flow and the thermal expansion water, 
consideration was given other probable factors which might not negligi- 
ble. These factors were grouped under Equation (1) and were con- 
sidered zero. 


Water budget evaporation 


The average monthly evaporation Lake Ontario for the period 1934 
1953 computed the water budget method was obtained substituting the 
values the normal change levels, runoff, inflow minus outflow, and pre- 
cipitation over the lake Equation (1). The values obtained are given 
Column Table above. Note that the minimum average monthly evapora- 
tion occurred June and the maximum average monthly evaporation oc- 
curred October. The reason for this will explained below. The total 
annual average evaporation for the period was 2.21 feet. 


Mass Transfer Method 


The velocity viscuous fluid the zone contact with solid boundary 
must have the same velocity the solid boundary. Thus, free stream 
air moving with the velocity past fixed solid boundary must have veloci- 
zero the boundary and velocity some distance away. Where 
velocity gradient found, the turbulent motion leads transport momen- 
tum across surfaces normal the velocity gradient. shearing stress 
exerted surface across which transfer momentum takes place. The 
relatively thin layer where the velocity gradient large called the bound- 
ary layer. Not only the character the fluid flow, laminar turbulent, 
important, but the nature the boundary has important function. The flow 
wind over Lake Ontario was considered turbulent and the surface 
aerodynamically rough. 

The general equation for shearing stress is: 


(2) 


where 
the shearing stress. 


the coefficient turbulent interchange momentum. 


the gradient horizontal wind velocity. 


The transfer water vapor through the atmosphere was considered oc- 
cur the same manner the transfer momentum. Thus: 


dz | 
where 
the evaporation. 


the coefficient turbulent interchange moisture. 


thi 


(2) 


(3) 


The specific humidity defined 


the gradient moisture concentration. 


(4) 
.378e 


where 
the vapor pressure the air 
the atmospheric pressure. 


Since the order one hundred times there little loss ac- 
curacy substituting 


(5) 


Since the coefficient turbulent interchange moisture was considered 
equal the coefficient turbulent interchange momentum, the 
evaporation was considered function the friction velocity, and 
the difference between the vapor pressure the air and the vapor pressure 
saturated air the water surface. 


where the friction velocity and equal /7/, 


Equations for the value the shear stress the wind over water under 
all conditions temperature stability have been determined.(3) using 
these equations, was assumed that when the water surface covered with 
waves different equations must used than when there are waves the 
surface. There appeared critical value for which the water 
surface becomes unstable, and waves form. 

Excellent data obtained from water-loss investigations Lake Hefner 
were used correlate measured evaporation with that computed Equation 
(6) above. The equation developed from the Lake Hefner data for evaporation 
inches per day which was considered applicable Lake Ontario was 


where 
centimeters/second. 


millibars and the vapor pressure the air measured 
meters elevation. 


The correlation factor Lake Hefner measured evaporation that com- 
puted the above formula was 0.90. The equation also appeared give good 
results for the Lake Mead studies. 

extremely difficult task obtain wind and temperature data for 
Lake Ontario. The friction velocity function the wind velocity over 
water given height and the atmospheric stability. The measure atmos- 
pheric stability used was the difference between the air and water surface 
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temperatures (TA TW). There are only two anemometers operating the 

year round near Lake Ontario where the wind velocities measured over land 

can adequately reduced over-lake wind velocities. Both these ane- 
mometer locations, Toronto and Cobourg, are Canada, and, unfortunately, 

the Cobourg anemometer was discontinued 1952. The recorded wind 

speeds were first reduced eight meter level wind speeds and then were in- 
creased over-water winds accordance with experience gained relating 
vessel weather observations Great Lakes vessels 

with land stations.(4) The factors which the over-land wind speeds were 
increased depended upon the atmospheric stability; these are given Table 


TABLE 


Ratio Over-Water Over-Land Wind Speeds, 


The vapor pressure differential (Ae) function the temperatures 
the air, the relative humidity the air, and the water surface temperature. 
The air temperatures and the dew points were derived arithmetically 


averaging the values recorded Toronto, Trenton, and Oswego, for the 
period 1937 through 1952. The water surface temperatures were derived 
from the work Freeman, Millar, and the Lake Survey, see Table 

The evaporation computed the mass transfer equation, Equation 
given Table 

Using the mass transfer method, the minimum average monthly evapora- 
tion occurred May and June and the maximum average monthly evaporation 
occurred September. The difference vapor pressures May and June 


1.50 1.21 1.00 
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TABLE 


Lake Ontario Water Surface Temperatures 


Degrees Fahrenheit 
Freem 


justed 


nan Temperatures 


September 
October 
November 


December 


the Great Lakes," report John Freeman. 


2usurface Temperatures the Great Lakes," Graham Millar, 
the Fisheries Research Board Canada, December 1952. 


Lake Survey, Detroit, Michigan 


shows that there was appreciable gradient moisture concentration and, 


consequently, net loss water mass. fact, some years, there was 


negative evaporation the months May and June. 


Discussion 


All data from known sources were used this study. However, the avail- 
able data were not the highest order, consequently, the computed evapora- 


tion subject error. New installations must established collect 
more complete and accurate meteorological and hydrological data the 


Month 
e 
January 
4 


TABLE 


Evaporation Feet Computed Mass Transfer Equation 


Over Water 
Meter Elevation Centimeters/second Millibars Foot/Month 
Centimeters/second 


January 
February 
March 
April 
May 

June 
July 
August 
September 
October 
November 


December 


Great Lakes area the most beneficial regulation plans are effected. 

comparison the evaporation computed the water budget and the 
mass transfer methods given Table 

The results compare rather favorably and the trend evaporation the 
same, although there are minor differences the amount net evaporation. 
The evaporation the spring months has been shown negligible because 
there appreciable gradient moisture concentration. The maximum 
evaporation occurs the fall. 

There definite evidence that the precipitation falling over the surface 
the Great Lakes considerably less than the precipitation falling the ad- 
jacent areas. The fact that the precipitation over Lake Ontario some 
per cent less than the precipitation recorded land stations throughout the 
basin means that the values evaporation are reduced the same amount 
when computed the water budget equation. result, the annual amount 
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TABLE 


Comparison Computed Evaporation, Feet 


Transfer 
Method 


Water Budget 
Method 


January 
February 

March 
April 

May 

June 

July 
August 
September 
October 
November 


December 


net evaporation from Lake Ontario computed both the water budget and 
the mass transfer methods the order 2.0 feet. 


Major Ira Hunt, Jr., Technique for Forecasting Observed Runoff from 
Large Areas, Lake Survey, (Detroit, 1958). 


Appointed the and Can. governments develop Great Lakes basic 
hydraulic and hydrologic data and establish the basis for future develop- 
ment such data these agencies, 


Capt. Ira Hunt, Jr., Effets vent sur les nappes liquides, (The Effects 
the Wind Liquid Surfaces), Univ. Grenoble, (Grenoble, France, 
1954). 


Major Ira Hunt, Jr., Relationships Between Over-Water and Over-Land 
Wind: Lake Erie, Lake Survey, (Detroit, 1958). 
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PROBLEMS CONCERNING USE LOW HEAD RADIAL GATES 


Thomas Rhone, A.M. ASCE 


ABSTRACT 


Three the principal hydraulic features low head radial gates are dis- 
cussed. These are (1) general discharge equation, (2) the effect gate seat 
location discharge capacity and pressure distribution along the spillway 
surface, and (3) side and bottom seals. 


INTRODUCTION 


The radial-type gate originated France about 100 years ago.(1) The 
earliest recorded use was the French Eagineer Poiree the Seine River 
1853. the 1860’s, another French Engineer, Mongel Bey, used cast- 
iron, radial-type gate the Delta Barrage the Rosetta Branch the Nile 
River. The Poiree gates were 28.7 feet wide 3.3 feet high; the Delta 
Barrage gates were 16.4 feet wide 16.7 feet high with the concave face 
turned toward the reservoir. 

The radial-type gate was patented the United States 1886. Since then, 
this type gate has become widely used and has increased size. Three 
114 feet wide feet high gates are used Horseshoe Dam Spillway near 
Phoenix, Arizona, and proposed use four 40-foot wide 52.5-foot high 
radial gates for reservoir control the new Glen Canyon Dam. 

Radial Tainter gates are probably the most widely used crest control 
gates. They are particularly well adapted crest control because their 
simplicity design, construction, and installation. The overall economy and 
efficiency this type gate result from the radial bearing. The thrust from 
the waterload carried two trunnion bearings where offers little resist- 
ance the operation the gate. This action can compared the hub 


Note: Discussion open until July 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1935 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. February, 1959. 


Hydraulic Research Engr., Division Eng. Lab., Bureau Reclamation, 
Denver, Colo. 
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wheel where all the force the rim transmitted the center. lift- 
ing the gate, there some force overcome the trunnion, and the side 
seals may offer considerable resistance the form friction. However, 
the opening radial gate requires less hoist capacity than slide gate with 
face bearing. addition, radial gates are more adaptable for automatic con- 
trol apparatus. 

Another desirable feature radial gate that needs gate slots. 
High head flow past indented slots has some installations produced cavita- 
tion damage the pier and spillway surfaces. One instance where gate slots 
were the apparent cause cavitation spillway and pier side walls 
shown Figure several other cases have been reported technical 
literature.(2) Preventive action has included the use gate followers fill 
the slots small gate openings, and the use offsets the pier side walls 
downstream from the slots. 

Despite the apparent simplicity and wide use radial gates, very little 
seems known about their operating characteristics. Specifically, there 
has been very little information published the discharge capacity, the ef- 
fect gate location pressure distribution the spillway face, downpull 
produced flow under the gate, and the type seals used, and other 
hydraulic problems. 

The purpose this paper discuss some the hydraulic features 
radial gates, including the side and bottom seals; present hydraulic data 
some representative gate installations; and create interest the remain- 
ing problems, with the hope that additional material will presented the 
readers this paper. Further, hoped that the tests described will 
stimulate hydraulic research radial gates those position 
hydraulic model prototype investigations. 

Material this paper has been limited the so-called low head radial 
gates used control flow over spillways. Top seal radial gates, where the 
normal operating head greatly exceeds the gate height, and the small gates 
used control canal flow, are not covered. 


Discharge Characteristics 


Discharge Determination 


Although the radial gate has been used for regulating flow for over 
hundred years, general method for predicting exact flow quantities has not 
been derived except for one type installation. When the gate used 
flat-bottom, rectangular section and the discharging jet supported the 
flat floor, the flow characteristics are known both from mathematical and 
experimental analyses. The discharge coefficients derived from these studies 
have been confirmed several experimenters and are accurate for their 
specific purpose.(3,4) 

Many specific discharge determinations are record; usually, when 
spillway has been model-tested, calibration curves for both free and gate- 
controlled flow have been obtained. From analysis these calibrations, 
the free flow discharge coefficients for almost any shape overflow section 
can determined.(5) However, general equation which will provide the 
exact discharges for radial-gate-controlled flow still unpublished. 

The lack general equation can probably explained the many vari- 
ables that affect the flow pattern. Some these are the spillway approach 
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FIGURE 


PARKER DAM SPILLWAY 


CAVITATION DAMAGE DOWNSTREAM 
FROM GATE SLOTS 


and include channel width and depth. Other variables include the shape the 
spillway profile; the geometry the gate, including the radius, width, height, 
trunnion location, and the location the gate seat; the location, spacing, and 
type piers; and the method determining gate opening and reservoir ele- 
vation. The combinations these variables are almost countless, and the 
effect any one them difficult determine from limited number 
investigations. 

brief review four methods which have been used determine the flow 
quantities under radial gates indicates the different approaches used at- 
tempts obtain general solution. should pointed out, however, that 
each provides only approximate discharges; accurate determination 
necessary, either model studies field calibrations should made. 


coefficient discharge listed function gate opening and reservoir 
head. The definition the symbols and suggested values for “C” are shown 
Figure From the graph, apparent that for values less than 
about 2.2 there wide range where “C” can vary much percent. 
The coefficient curve shown the heavy line was derived from hydraulic 
model studies the spillways listed Figure The different crest shapes 

and gate arrangements probably account for some the spread. However, 
known that where the head orifice small comparison the 
height the orifice, there appreciable difference between the discharge 
obtained using average head and the discharge obtained taking into 
consideration the variation head. Depending the actual values the 
orifice height and head, the differences discharge may large small. 
Thus, the region below 2.2 transition region, region notoriously 
difficult evaluate all hydraulic flow problems. 

second method determining discharge quantities takes into considera- 
tion the angle lip gate horizontal and the angle described the 
intersection horizontal line through the pivotal point the gate and the 
radius the gate drawn from the same point the center the gate when 
closed position, Figures The value the coefficient dis- 
charge for controlled flow related the known coefficient discharge for 
free flow for different values and ratios gate openings depth 
flow hy. Plots for two different values are shown Figures and 
The rate flow then obtained applying the formula 


3/2 
Qe hz -ho 


where vertical distance from crest water surface 


= COEFFICIENT OF DISCHARGE 


vertical distance from bottom gate water surface 
width gate 


The coefficients were derived from the results many model studies. 
The third method similar the second that relates the controlled 
flow discharge coefficient the free flow coefficient, Figure The effect 
gate radius, trunnion location, and gate opening also recognized having 
effect the coefficient and evaluated the parameter The 
s | 
values for “C” were obtained from analyses calibrations model spillways. 
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DISCHARGE UNDER RADIAL GATE 
FIGURE 5 


The discharge formula for this method is: 


Values for and 1/sin are given Figure 

The fourth method, Figure takes into consideration the crest shape, lo- 
cation gate seat, the angle formed the tangent the gate lip and the 
tangent the crest curve the nearest point the crest curve, and the net 
gate opening which the shortest distance from the gate lip the crest 
curve. From the geometry the gate the angle determined and from the 
curves Figure coefficient discharge determined. These values 
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for “C” were also obtained from model investigations. The discharge then 
computed from the formula: 


where discharge cfs 
net gate opening 
crest width 
head the center the opening 


Table shows Canyon Ferry Dam Spillway discharges computed from the 
four methods described, compared with the discharges obtained from model 
and prototype calibration. The Canyon Ferry Dam Spillway radial-gate- 
controlled spillway located the Missouri River about miles east 
Helena, Montana. The spillway was calibrated model studies 1953, and 
field discharge measurements were obtained 1956.(6) 

The first and fourth methods give percent higher discharges than 
actually measured. The second and third methods give discharges from 
percent over 7.5 percent under the measured values. The small gate 
opening shows greater variation than the larger openings. 

The obvious conclusion that can drawn from Table that satisfac- 
tory method for computing the discharge through partially open radial gates 
has not yet been derived. method computing discharges required which 
will better agreement with model tests and field calibrations. Examina- 
tion the four discharge equations shows that three them the discharge 
varies and the fourth more nearly correct method 
might vary this exponent with changes The change from free 
weir flow orifice flow, vice versa, would, thereby, taken into ac- 
count. The transition range occurs when less than 2.2, Figure 
When the ratio greater than 2.2, true orifice flow occurs. Thus, the 
free flow range the exponent for the head the discharge formula three- 
halves; the transition range, the exponent should vary between three-halves 
and one-half; for orifice flow, the exponent one-half. 


Effect Gate Position Discharge Coefficient 


1953, hydraulic model investigations were made determine the effect 
pressure distribution along the spillway profile caused changing the lo- 
cation the gate seat.(7) 

For this study, crests were investigated. With each crest, gate 
seat locations were tested. The locations were (1) the crest, (2) and (3) 
inches vertically below the crest both upstream and downstream from the 


openings and heads are measured vertically above gate 


seat. 
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crest axis, (4) inches below the crest the downstream side, and (5) the 
design gate seat location when did not correspond one the other four. 

The five crests and the model scale ratio which each was studied are 
(1) Ross Dam 1:24.14, (2) Altus Dam 1:20.70, (3) American Falls Dam 1:18.35, 
(4) Bhakra Dam 1:50.56, and (5) Canyon Ferry Dam 1:42.92. The crest pro- 
files are shown Figures 12. 

part these tests, data necessary compute the coefficient dis- 
charge for several gate openings with each gate location were obtained. The 
coefficients were computed using the first method described. was found 
that when the gate seats either upstream from the crest the crest, the 
coefficient discharge higher than when the gate seats downstream from 
the crest. This was true for all crests; some cases, the difference was 

Some the difference might explained the fact that the head and gate 
opening were measured from the gate seat instead the crest. However, the 
head and gate opening values are the same whether the gate seats inches be- 
low the crest upstream downstream direction from the crest. all 
cases, the coefficients for the upstream position were always higher. ex- 
ample these coefficients shown Figure 


Conclusions and Recommendations 


Four different discharge equations have been used compute the flow 
quantities under the radial gates the Canyon Ferry Dam spillway. The 
discharge coefficients used each equation were based different model 
studies; each equation placed emphasis different geometric features the 
structure. one equation, the head was raised the one-half power; the 
others, the head was raised the three-halves power. The flow quantities 
obtained from these methods were compared with actual model and prototype 
measurements; the computed flows were found vary from the measured 

All the methods described would probably give discharges sufficiently 
accurate for design purposes for discharge determinations where water 
has little value and rough estimate the flow would suffice. However, when 
necessary have accurate discharge determination for flood routing 
measure valuable irrigation water, the western United States, 
either model prototype calibration necessary. 

accurate formula could derived combining theory and data 
sis from coordinated hydraulic investigations model prototypes. The 
amount existing data available from hydraulic laboratories might provide 
sufficient material for preliminary analysis and these data were supple- 
mented models designed for discharge studies, the effect 
each variable could probably evaluated. Such models would have 
sufficiently large that changes the geometry the gates, piers, and over- 
fall section would make measurable changes the flow. The crest length 
should sufficiently long allow installation several piers, and provision 
should made test several crest shapes, pier shapes, and gate sizes. 
would desirable incorporate adjustable weir the model headbay, 
which would maintain constant reservoir level facilitate the flow measure- 
ment the transition range between free flow and orifice flow. Standard 
measuring devices for determining gate openings, discharge quantities, and 
heads the crest would adequate. With such model, the effect any 
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the variables such approach depth, pier shape, trunnion location could 
evaluated. 


Pressure Characteristics 


standard spillway crest usually designed conform the shape 
the lower nappe downstream from sharp-crested weir. This profile 
called the datum shape.(5) When gates are used the spillway, the 
same crest shape used from the upstream face the gate seat. Down- 
stream from the gate seat, the crest profile sometimes made corres- 
pond the theoretical trajectory jet flowing from 1-foot high rectangu- 
lar orifice the design head. For the same design head, the crest designed 
for free flow will steeper than crest designed for 1-foot gate opening. 

Several factors other than hydraulic phenomena are usually considered 
when determining the location the gate seat. the gate seats the crest, 
the theoretical crest curve will comparatively flat and the spillway will 
contain considerably more concrete than usually needed for stability. 
the gate seats below the crest, the profile will steeper, but higher gate 
necessary retain the same storage pool. Structural and architectural 
considerations sometimes dictate the gate seat location. bridge spans 
the spillway, might necessary place the gate downstream from the 
crest provide space for the gate lifting mechanism provide clearance 
when the gate raised. The gate might have located downstream in- 
sure that the trunnion will above the water surface. 

Many hydraulic designers believe that when the gate seated the crest 
there tendency for the jet spring free the crest high reservoir 
elevations and small gate openings; this may result subatmospheric pres- 
sures and possible cavitation damage the spillway. The designers also be- 
lieve that placing the gate seat downstream from the crest results the flow 
being directed downward; the tendency the jet spring free the crest 
lessened. However, observations prototype spillways indicate that cavita- 
tion damage the spillway surface occurs downstream from the region where 
the control gates have any effect the flow. damage spillway crests 
the vicinity the gates has been found which can attributed the crest 
shape and gate seat location. 

realistic solution these problems and conflicting requirements 
needed. The following discussion may help clarify the problems. 


Effect Gate Location Pressure Distribution Along 
Spillway Profile 


The laboratory investigations used determine the effect gate location 
discharge coefficients were also used determine the effect the gate 
position the pressure distribution along the spillway surface.(7) For the 
pressure tests, piezometers were installed along the crest profile. One 
piezometer was placed each tested gate seat and the others were equally 
spaced downstream along the spillway profile. The five crests used the 
investigation provided wide range steepness and gate seat locations. 

For convenience, each the crests will discussed separately giving 
description the crest and the results the tests. The crest profiles, 
datum shapes, piezometer locations, and pressure profiles for representa- 
tive gate openings are shown Figures 10, 11, and 12. 
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Ross Dam, Figure 

The design gate seat inches below the crest. The datum shape 
slightly steeper than the design shape. comparison the design shape and 
datum shape indicates that the overfall section from the gate seat downstream 
was designed from the theoretical trajectory jet from orifice. From 
the reasoning previously outlined, pressures along the profile should increase 
when the gate seat downstream from the design location and reduce when 
the gate seat upstream. 

The pressure profiles Figure show that for the same gate openings 
and heads, the pressures for all gate seat locations downstream from the de- 
sign gate seat are practically the same. However, when the gate seated 
upstream from the crest, the piezometers from the crest down the second 
piezometer below the design gate seat showed pressures lower than were ob- 
tained for the other gate seat locations. Another trend which may ob- 
served from the pressure profiles that the zone subatmospheric pres- 
sures moved upstream the gate seat moved upstream. However, for all 


gate positions, the downstream piezometer (No. showed pressures near 
atmospheric. 


Altus Dam, Figure 

The design gate seat the crest. The datum shape the same the 
design shape, indicating that the profile was obtained from the lower nappe 
shape sheet water passing over sharp-crested weir. According 
reasoning, the pressures along the profile should have the highest values 
when the gate seat inches below the crest and should become smaller 
the seat moved upstream. The pressure profiles Figure show this 
true that for both downstream gate seat locations the piezometers 
downstream from the gate seat indicate nearly atmospheric pressures for all 
gate openings. Also, when the gate seats upstream from the crest, the piezo- 
meters from the gate seat Piezometer No. indicate subatmospheric pres- 
sures. The piezometers downstream from No. all show near atmospheric 
pressures. Generally speaking, Piezometers No. through showed nearly 
atmospheric pressures for all gate seat locations all gate openings. Piezo- 
meters No. through were affected the gate seat location and the pres- 
sure was reduced each piezometer the gate seat was moved upstream. 


American Falls Dam, Figure 

The design gate seat the crest. The datum shape considerably 
steeper than the design shape, indicating that the profile from the crest down- 
stream was obtained from the theoretical trajectory jet flowing from 
orifice. The pressures should near atmospheric for small gate openings 
when the gate the design location, for downstream gate seat locations 
the pressures should higher, and for upstream locations the pressure 
should lower. The pressure profiles Figure indicate that perform- 
ance agrees with the reasoning. The only pressures significantly subatmos- 


pheric were those the crest piezometer when the gate seat was upstream 
from the crest. 


Bhakra Dam, Figure 

The design gate seat 2.43 feet below the crest the downstream side. 
The datum shape steeper than design shape. For this crest, the highest 
pressures should occur with the gate the design location and they should 
become successively lower the seat location moved upstream. Also, 
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the piezometers should indicate near atmospheric pressures when the gate 
seats the crest. The pressure profiles Figure show this es- 
sentially true. The piezometers downstream from the gate seat show slight 
reductions pressure the seat moved upstream and when the gate seats 


the crest and upstream from the crest, the pressures are little above and 
below atmospheric. 


Canyon Ferry Dam, Figure 

The design gate seat inches below the crest. The datum shape 
steeper than the design shape. The relative shape the two profiles indi- 
cates that the profile downstream from the design gate seat was derived from 
the theoretical trajectory jet flowing from orifice. Therefore, the 
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PIEZOMETER ond gate openings measured vertically above crest 
(2) GATE SEAT AT PIEZOMETER 4, 12" BELOW CREST ~~---==-Gate open 143% of design head 
————— Gate open 426% of design head 
Gote open 572% of design head 
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pressures downstream from the gate should near atmospheric when the 
crest the design location, higher when the gate seat downstream and 

lower the gate seat moved upstream. The pressure profiles Figure 
indicate that the pressures increase when the gate seat downstream 
from the design seat, and decrease when the gate seat upstream. However, 
for gate openings 8.8 and 23.5 percent the design head, the pressures 
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downstream from gate seat are practically the same for any gate seat 
location. 


Conclusions 


The pressure tests show that the gate seat location has minor effect 
the profile pressures when the gate seat downstream from the crest. 
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When the gate seats upstream from the crest, there sharp reduction 
the pressure the crest, particularly the small gate openings. 

The theory that gate seat downstream from the crest will give down- 
ward direction the flow passing beneath the gate seems confirmed 
these tests. Although the increase pressure was small, was noticeable 
for all crests. 

The lowest pressures occurred the two crests that most nearly corres- 
ponded with the datum shapes, Figures and and generally occurred for 
gate opening about percent the design head. For the other three 
crests, Figures 10, 11, and 12, the lowest pressures occurred the 
gate openings. The lowest pressure recorded was about percent the de- 
sign head below atmospheric. 


Gate Seals 


Radial gates are usually provided with side and bottom seals. some 
installations the side seal bears against steel plate embedded the face 
the pier; this provides smooth surface for the contact and prevents exces- 
sive wear the sealing strip. The bottom seal contacts metal plate set 
flush with the spillway surface; shallow depression about 3/32 inches 
the metal sill plate acts seating surface. 

1952, limited number tests were performed the Bureau 
Reclamation Hydraulic Laboratory determine the effectiveness the 
gates seals being used.(8) Because limitations the test rig, sliding 
the side seal the wall plate corresponding the gate opening closing 
was not possible. The seal could moved direction normal the wall 
plate, determine the force required seat the seal. When the seal seated, 
leakages were too small conveniently measured, visual 
the leakage sufficed for the tests. 

Two tests were made with each seal. The first test was determine the 
lateral thrust per lineal foot seal. the lateral thrust was within reasona- 
ble limits, the seal was tested the hydraulic rig determine its action un- 
der pressure. 

The lateral thrust tests were made with wooden mock-up the clamp 
plates radial gate, Figure 13. The thrust was determined lifting the 
seal and off the laboratory floor with spring scale and steel bar placed 
under and perpendicular the seal axis. The seal was lifted point mid- 
way between two pieces paper placed inches apart the floor under the 
seal. When the paper could slipped free, the thrust load the floor was 
considered relieved and the spring scale load recorded. The frictional drag 
was computed from this thrust value. 


Belt-type seal 

One proposed side seal was made from fabric-reinforced rubber belting. 
seal inches wide and feet long was assembled the wooden mock-up 
for lateral thrust tests. The lateral thrust was determined for various seal 
positions. Zero seal position occurs when the space between the end the 
gate and the wall plate three-quarters inch. The tests showed that 
the zero seal position, the lateral thrust was pounds per lineal foot seal. 
This thrust was high; and with the addition waterload 40-foot head, 
the frictional drag would pounds per inch seal. view the high 
drag, hydraulic tests were performed. 
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(A) General view the 3/8 inch foot long 
belting seal bolted 40-foot radius. 
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(B) Close-up view the above assembly. 


FIGURE 


RADIAL GATE SEAL TESTS 


Wooden Mock-up Rubber Belting Side Seal without Brass Shoes 


al. 


Belt-type seal with brass shoes 

was decided investigate the same seal with brass shoes bolted the 
belting reduce the frictional drag, Figure 14A. The lateral thrust zero 
seal position for this seal was 158 pounds per lineal foot; but because the 
smaller coefficient friction with the brass shoe, the frictional resistance 
40-foot head was pounds per inch seal. 

Hydraulic tests showed that when the seal was moved away from the wall 
plate, the seal shoe slid downstream and remained flat contact with the wall 
plate. However, when the gate was moved nearer the wall plate, the brass 
shoe did not slide upstream the wall plate, and the rubber belting buckled. 
Since this action was independent the reservoir head, the seal was con- 
sidered unacceptable. 


Molded-type belt seal 

means reducing the lateral thrust the belt type seal, was 
proposed that the belt seal molded angular shape, Figure 14B. The 
lateral thrust this seal was negligible zero seal gap. The frictional drag 
the wall plate was computed about pounds per lineal inch. 

Hydraulic tests indicated that the brass shoe did not slip upstream freely 
the seal gap was decreased. The seal buckled when the gap was decreased 


1/4-inch, and the sealing lip raised slightly from the seat. was concluded 
that the seal operation was only fair. 


Angle seals 

Two types angle seals were investigated. One was all-rubber seal, 
and the second was similar, except that was fabric-reinforced, Figure 15. 
The lateral thrust with these seals would depend how they were mounted 
the gate sides. Because this, tests were performed determine 
what head these seals would seal and unseal for various seal-to-side wall 
spacing. addition, identical tests were made for the all-rubber seal with 
brass bar the sealing edge. 

For the all-rubber angle seal, without the brass bar, the water head re- 
quired for sealing 1/2-inch gap was feet, while unsealing occurred 
8-foot head, Figure 16. With the brass bar, these values were and 
feet. The fabric-reinforced seal sealed 19-foot head and unsealed 
6-foot head. The computed frictional drag for these seals was about 
pounds per lineal inch without brass bar and 6.5 pounds per lineal inch with 
brass bar. 

Although the angle seal with the brass bar appeared have better sealing 
and unsealing characteristics, well lower frictional drag, allowed 
slightly greater leakage than the all-rubber seal and probably would less 
desirable for this reason. 

Two seals were made that included rubber-to-metal sealing surfaces 
brass member reduce frictional drag, Figure 17. The seal with the rub- 
ber flaps (Figure 17A) was generally unsatisfactory, since the head required 
obtain water seal normal gaps was higher than normally found for 
spillway gates. 

The second seal, shown Figure 17B, incorporated more stabilized 
sealing lip. 1/2-inch gap, this seal sealed 20-foot head and un- 
sealed 4-foot head, Figure 18. The frictional drag with this type seal 
40-foot head and zero gap approximately pounds per lineal inch. 
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(A) View of the rubber belting seal with brass shoes. 


(B) View of the belt type, molded angle seai, with brass shoe, 
cut from a fabric reinforced angle seal shown on Figure $a. 
FIGURE 
RADIAL GATE SEAL TESTS 
Belt Type Side Seals with Metal Shoes 


‘ed 
| 


February, 1959 


(A) View the angle seal with fabric reinforcement. 


(B) View the all-rubter angle seal with brass bar added 
reduce frictional drag wall plate. 
FIGURE 
RADIAL GATE SEAL TESTS 


Views the Molded Angle Side Seal with Fabric 
Reinforcement and Brass Bar 
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sad 
(A) View the angle seal with brass bar plus rubber sealing flap. 
" (B) View of the angle seal with rubber removed from the lower leg & 
and brass plate added. 
FIGURE 
RADIAL GATE SEAL TESTS 
Views the Angle Side Seal with Brass Bar plus 
Rubber Sealing Flap and with Brass Plate. 
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Conclusions 


Based the results these tests and the computations for frictional 
drag, side seal design should include the features listed below. 
design incorporating these features shown Figure 19. 


The sealing force should not dependent the water head. 
rubber-to-metal seal preferred obtain the least water leakage. 
The portion the waterload carried the rubber seal should low 
reduce the frictional drag. 
The heel end the seal should mounted that will not contact 
the wall plate and increase the frictional drag. 
The major part the seal waterload should carried through metal 
member the wall plates minimize the frictional drag the seal. 


Radial Gate Bottom Seal Tests 


The hydraulic test rig used the side seal tests was modified simulate 
bottom seal arrangement, Figure 20. The seal seat was made with 
depression 3/32-inch deep and inches long measured parallel the long 
dimension the seal represent damaged seating surface which might 
encountered the field. The angle the seal seat relative the seal was 
similar that field gate. Two all-rubber seal samples were available 
for test; each one was rectangular section. One was made stock with 
Type Shore durometer hardness 69; the other was made stock with 
hardness 38. 
Both seal samples were tested the hydraulic rig heads 195 feet 
‘of water. This test consisted compressing the seal against the simulated 
gate seat and increasing the hydraulic head 195 feet water. Both seals 
appeared seal satisfactorily, since there was leakage the depression 
the gate seat. However, heads feet more, the soft rubber seal 
had tendency creep under the downstream clamp plate. 

The principal test conducted these two seal samples was determine 
the force required compress the seal. These data were desired for gate 
design purposes assure sealing depressions the seal seat. This in- 
formation was obtained Universal-type testing machine determining 
the vertical force required compress 1-foot wide section the seal 
against wet and dry seats. The wet seal data were obtained submerging 
the simulated seal seat pan water. plot the compression force, 

pounds per lineal foot seal for wet and dry seat, shown for both seal 
samples Figure 20. 


Conclusions 


The rectangular-type radial gate bottom seal tested was satisfactory. The 
softer rubber material would probably form more effective seal against 
surface irregularities the seal seat. The tests showed that for surface 
irregularities approximately 3/32 inch deep, the soft material would form 
seal applied force pounds per lineal foot seal; the harder ma- 
terial would require force 140 pounds per lineal foot. 
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NOTES 
water load across the 
mension carried the gate, 
the remainder water load 
carried the wall plate. 
2.Brass pipe may slightly curved 
during assembly fit gate 


tal 3.At higher heads,seal leg will 

tend curve,and thus releave 

the frictional drag rubber 

sealing lip wall plate. 

4.Drawing full scale. 
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Holes, Centers 
seal wall plate. 

RADIAL GATE SEAL TESTS 
PROPOSED SIDE SEAL WITH BRASS PIPE 
AND RUBBER SEALING LIP 
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COMPRESSION PER FOOT SEAL 
FORCE REQUIRED TO COMPRESS RADIAL GATE BOTTOM SEAL 
RADIAL GATE SEAL TESTS 
ARRANGEMENT FOR HYDRAULIC AND COMPRESSION TESTS BOTTOM SEAL 
FIGURE 
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BOUNDARY-SHEAR STRESS UNSTEADY 
TURBULENT PIPE FLOW 


Carstens,! A.M. ASCE and John Roller,2 J.M. ASCE 


INTRODUCTION 


There are two reasons for the preponderance experimental and theoreti- 

cal knowledge concerning steady flow circular pipes compared with un- 

steady flow. First, problems involving steady quasi-steady flow pipes 

arise repeatedly the majority engineering works. Second, experimental 

research unsteady flow has been delayed until the development reliable 

dynamic-pressure-measuring instruments. 
this paper, theoretical analysis boundary-shear stress for unsteady 

turbulent flow pipe given, followed report experiments with un- 

steady flow smooth pipe. 

There good reason restrict the study turbulent flow. laminar 

flow, the instantaneous mean velocity and instantaneous velocity profile (hence, 

the boundary-shear stress) are dependent upon the history the motivating 

force piezometric-pressure gradient. result this historical de- 

pendence, each case unsteady laminar flow unique. Conversely, tur- 

bulent flow, the strong lateral diffusion the turbulent eddies tends elimi- 

nate the dependence the velocity distribution upon the history the motion. 

this hypothesis justified, then the dependent flow variable, boundary- 

shear stress, function instantaneous values the flow variables, 

mean velocity and acceleration (or piezometric-pressure gradient). Thus, 

greater degree generalization results appears possible for un- 

steady turbulent flow pipe than for laminar flow. 


Note: Discussion open until July 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1945 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. February, 1959. 
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NOMENCLATURE 


cross-sectional area pipe; 


inside diameter pipe; 

Darcy-Weisbach coefficient boundary-drag coefficient; 

acceleration gravity; 

piezometric head; 

piezometric head the reservoir; 

length pipe 

reciprocal the exponent the power law velocity distribution; 
piezometric pressure; 

discharge volume rate flow; 

radial coordinate measured from the centerline the pipe; 
pipe radius; 

fluid density; 

shear stress; 

boundary-shear stress; 


time, measured from the instant initiation the flow; 


velocity the pipe; 

mean velocity the pipe; 

“observed” velocity determined from the jet coordinates; 
volume fluid element within the pipe; 


axial coordinate measured from the upstream end the pipe; 


subscript designating unsteady flow; and 


subscript designating steady flow. 


Theoretical Analysis 


The linear-momentum equations are the beginning point rational analy- 
sis boundary-shear stress. uniform unsteady flow pipe, the velocity, 
everywhere axial and function and only. Similarly, the shear 


acceleration, the piezometric pressure, function only and 


The fluid element, length dx, extends the pipe wall, shown 


Figure The linear-momentum equation for this element 
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Figure Fluid Element. 


The terms the left side the equation are the external forces the 
fluid element. The first term the combined pressure and weight force 
the direction. The second term the boundary-shear force. 

The first term the right side the equation the difference momen- 
tum flux through the end areas the fluid element given time. Since the 
equation being applied uniform flow, this term zero. The existence 
turbulent velocity fluctuations has effect the magnitude this term. 

The second term the right side the equation the time rate change 
linear momentum the particles within the interior the element. Since 
plifications are possible, 


Again the existence turbulent velocity fluctuations has effect the 
magnitude the term. 


The simplified linear-momentum equation for the fluid element Figure 


(1) 


linear-momentum equation can derived similar manner 


with annular cylindrical element radius and thickness dr. The result- 
ing linear-momentum equation 
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Equation (2) function the velocity distribution. contrast, the cor- 
responding term for steady flow constant regardless the velocity 
distribution. 
the magnitude the acceleration small, reasonable assumption 
that the velocity distribution the same for the unsteady flow for the cor- 
responding steady flow. The basis for the assumption the strong lateral 
diffusive effect the turbulent eddies. The process can visualized 
considering the chronological sequence events initially steady flowis 
accelerated another steady-flow state. Any deviations from the steady- 
state velocity distribution which occur during acceleration will rapidly 
removed the turbulent-eddy diffusion the final steady state. addition, 
the turbulent-eddy diffusion also present during the acceleration. Thus, 
for turbulent flows undergoing moderate acceleration, equality the 
state and steady-state velocity distributions rational assumption. 
The power law velocity distribution for turbulent flow smooth pipe 
used the subsequent analysis. The power law velocity distribution 


(1). 


The value seven for pipe Reynolds numbers less than The 
value increases with increasing Reynolds number, becoming ata 
Reynolds number 3.2(106). 

Equation (2) integrable with respect upon substitution the power 
law velocity distribution. The piezometric-pressure gradient eliminated 
means Equation (1). The resulting expression for shear-stress 


and for shear-stress gradient 


Values and are given the Table for The 
The unsteady-flow shear-stress Equations (3) and (4) would more 
meaningful related the well-established expressions for shear stress 
steady flow. This relation possible virtue the equality the shear 
stress and the shear-stress gradient the centerline for both steady and un- 


steady flow. Since similar velocity distributions were assumed steady and 
unsteady flow, the deviation unsteady-flow shear-stress distribution from 


linear 
accele 
stress 
the 

betwee 
distril 
axis, 
infi 
Thus, 


But 


Introd 


Simila 


Z 
add: 
steady 
Values 
embod 
Fig 
ducing 
The 


linear must solely the result acceleration. all fluid particles were 
accelerated equally, rigid-body motion, the internal stresses (shear 
stresses) would not influenced the acceleration. Therefore, the change 
the shear-stress distribution must the result relative acceleration 
between adjacent particles the cross section. Inasmuch the velocity 
distribution symmetrical about the pipe axis and continuous through the 
axis, the relative acceleration between particle the axis and particle 
infinitesimal distance from the axis will also infinitesimally small. 
Thus, consequence the assumption identical velocity distributions 
unsteady and steady flow, the shear-stress gradients the axis are equal, 


addition, the shear stress zero the pipe axis for either steady un- 


But steady flow 
Introducing these concepts into Equation (4), 
Similarly, for Equation (3), 


Values are given the Table for The physical concepts 
embodied Equations (4), (5), and (6) are more easily visualized referring 
Equation (5) the desired relationship between the boundary-shear stress 
unsteady flow and that steady flow. Rearranging Equation (5) and intro- 
ducing the boundary-drag coefficient, results 


The boundary-drag coefficient, influenced two ways the magni- 
ind tude the Reynolds number. First, function only the Reynolds 
m | 


for flow smooth pipe. Second, the value which characterizes 
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Figure Shear Relationships Unsteady Turbulent Flow. 


the shape the velocity profile, also function the Reynolds number. 
The numerical value, 0.449, Equation (7) for This numerical value 
0.391, 0.346, and 0.310 for eight, nine, and ten, respectively. 

Thus, the theoretical analysis indicative that the boundary-shear para- 
meter, unsteady and turbulent flow smooth pipe independent 
the history the motion and dependent only upon the corresponding steady- 
state value and the instantaneous value the acceleration parameter, 
(D) V2. Equation (7), which mathematical expression for this 


dependence, was formulated using the hypothesis that the velocity distribution 
similar for both steady and unsteady flow. For moderate acceleration, the 
hypothesis appears valid, but experimental results are desirable 


verify and establish the limits the hypothesis. 
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TABLE: and for 


Equation Equation (4) Equation (6) 
-0.1123 
-0.011 -0.0962 0.0002 
0.2 -0.0200 -0.0861 0.0025 
0.3 -0.0279 -0.0706 0.0058 
-0.0340 -0.0534 0.0109 
0.5 -0.0385 -0.0320 0.0176 
0.6 -0.0403 0.0261 
-0.0396 0.0256 0.0378 
0.8 -0.0349 0.0706 0.0549 
0.9 -0.0224 0.1414 0.0786 
1.0 1.0000 0.1123 


Experimental Study 


Whether steady unsteady flow, experimental determination boundary- 
shear stress indirect. unsteady flow, this determination, means 
Equation (1), complicated the necessity measuring p*(x,t) and V(t). 

these experiments, the function V(t) was determined analysis 
motion-picture record the free jet the pipe outlet. Commercially avail- 
able pressure transducers were used determine p*(x,t). 

The boundary and initial conditions the experiment were follows: 

(a) the smooth, straight, circular pipe was attached large reservoir 
water which the piezometric head was maintained constant; and (b) the 
flow was initiated the instantaneous opening downstream valve. Since 
the fluid was rest prior the opening the valve, the initial flow was 
laminar throughout the pipe. With the passage time, spots turbulence 
were generated. These turbulent spots increased size with increasing 
time until the flow was turbulent throughout the pipe. 


Experimental Equipment and Technique 


The overall arrangement the equipment shown Figure The test 
section was extruded 1/2-inch brass pipe. The experimentally determined 
piezometric-head gradient six steady-flow runs was within the range, +0.1 
per cent -1.08 per cent, the piezometric-head gradient computed from 
the Blasius equation. This close agreement means that the pipe was smooth 
and that the diameter, measured means micrometer, was correct. 
Piezometers were spaced intervals diameters along the pipe. The 
signals were amplified electronically and were recorded 
Sanborn two-channel oscillograph. The configuration the downstream 
jet from the horizontal pipe was recorded means motion- 
picture camera. 
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Experimental Data 


PIPE FLOW 


The oscillographic records pressure were replotted instantaneous 
piezometric-head lines order delineate the and turbulent-flow 
regions. The succession piezometric-head lines typical run are shown 
Figure The abscissa Figure the distance along the pipe axis 
from the reservoir. The ordinate The position the zero value 
the ordinate lowered with increasing values order 


avoid the superposition many piezometric-head lines. The dashed line 
Figure delineates the regions laminar and turbulent flow. Laminar-flow 
regions are above the dashed line and turbulent-flow regions are below the 
dashed line. 

With the free-fall equations, the “observed” velocity-time relationship 
the downstream jet was determined from the coordinates and from the time 
recorded the film strip. Two assumptions are implicit the use 
the free-fall equations. First, air resistance negligible the jet. 
Second, the jet composed individual masses which each fluid particle 
has the same linear momentum the point measurement. Two experimen- 
tal observations substantiated the validity the assumptions, follows: 

(a) the “observed” velocity-time relationships were identical whether the 
point measurement was 0.5 below the pipe outlet whether the point 
measurement was 2.50 below the pipe outlet; and (b) all the film strips, 

discrete masses fluid appeared forming the jet the point 
measurement. However, the “observed” velocity, V', not the mean veloci- 
the pipe. 

The mean velocity, the pipe less than the “observed” velocity, V', 
determined from the jet coordinates. The difference the two velocities 
due the radial variation velocity within the pipe. The radial variation 
velocity the pipe outlet the result the boundary-shear force the 
fluid passes through the pipe. the fluid passes the pipe outlet, the boundary- 
shear decreased suddenly negligible value and the internal shear 
forces within the jet tend eliminate the non-uniform velocity distribution. 
The linear momentum the jet greater the pipe outlet than the veloci- 
distribution were uniform. The mass fluid passing the pipe outlet 
time PQat. The linear momentum this fluid mass PQVat. 
The linear-momentum-correction coefficient, measure the non- 
uniformity the velocity distribution the fluid passes the pipe outlet. 
the absence any horizontal external forces the linear momentum this 
mass constant the jet. the point measurement jet coordinates, 
the linear momentum Therefore, the desired quantity 
equal Only the value the “observed” velocity, V', can de- 
termined from the film strip the jet. Other methods had employed 
determine C,, 
The experimentally determined values laminar flow were in- 
corporated into the solution the Navier-Stokes equations order deter- 
mine the velocity distribution the pipe outlet function time. 
noted Figure that the flow was laminar the pipe outlet for 
<1.38 even though turbulent flow existed portions the pipe 
much earlier time. The heat diffusion equations analogous the Navier- 
Stokes equations have been solved with analogous boundary and initial condi- 
tions and with the analogous piezometric-head gradient arbitrary func- 
tion time.(2) The piezometric-head gradients were determined from 
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Figure and were replotted Figure simple analytical function, 

was selected being good approximation for the experi- 
mentally determined piezometric-head gradient laminar flow. course, 
had determined empirically for each the eight runs. This function 
was then incorporated into the solution the Navier-Stokes equations for the 
determination C,, function time. Because the length the 
laminar-flow region was less than the distance between pressure measure- 
ment stations when the piezometric-head gradient 
laminar flow could not experimentally determined from Figure when 

Vgho/L Thus during the interval 1.00 1.38, had 
determined some other manner. 

Figure was utilized obtain the value Cy, this time. The discon- 
tinuity occurs the time the turbulent flow first passes the outlet. The dif- 
ference velocity distribution between the turbulent flow and preceding 
laminar flow results discontinuity linear momentum the laminar 
and turbulent jets this time even though continuous. The value 
the turbulent flow was assumed 1.02 based upon the seventh-root 
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equal for both laminar and turbulent flow this time, the experimental 
point, shown circle Figure was readily determined. The dashed 
line Figure interpolation between the two methods. 
The resulting V-t function shown Figure 6.* 


Experimental Results 


The boundary-shear stress was determined indirectly means 
Equation (1) utilizing the data shown Figures The results are 
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Figure "Observed" Velocity, Mean Velocity, and Linear-Momentum 
Coefficient Typical Run. 


more complete discussion the method analyzing the data contained 
the paper “Transition from Laminar Turbulent Flow Pipe,” 


M.R. Carstens, Journal the Hydraulics Division, Proc. A.S.C.E., Paper 
1450, December, 1957. 
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presented Figure the form suggested Equation (7). Results ex- 
periments conducted Massachusetts Institute Technology(3) 
unsteady-flow water tunnel are also shown Figure 


% 
Analysis Results 


There are two obvious explanations for the scatter the experimental 

results shown Figure 

First, the scatter the experimentally determined values indication 
the precision determining the drag coefficient, The p-t function 

could measured within probable error per cent regard pres- 

sure. view the complexities determining the V-t function from the 

motion-picture data, probable error per cent would seem 
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Figure Experimentally Determined Values Boundary-Shear 
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optimistic. further source possible error the graphical determina- 
tion dV/dt from graphs similar Figure Since none these possible 
errors compensating, the entire scatter the experimental results might 
attributed logically experimental error. 
Second, the unsteady-flow drag coefficient, may function the 
deduced from Equation (7). However, the basic assumption 
the derivation Equation (7) that the history the turbulent flow does 
not influence the instantaneous flow variables. Thus, history the 
motion important, the scatter the experimental results Figure might 
the result omitting the parameters necessary define this history. 
this regard, Daily, al(3) noted historical effect during deceleration 


not during acceleration. 


The experimental results are not sufficiently precise either prove 
disprove the validity the assumption that the history turbulent flow does 
not influence the instantaneous flow variables. The derivation, leading 
Equation (7) for the drag coefficient unsteady flow considered ra- 
tional. 

Both the results this study and those Daily, al,(3) indicate that the 
current engineering practice assuming that =f, acceptable for turbu- 
lent flow pipes. the first place, the relative importance the boundary- 
drag force decreases unsteady flow compared steady flow while the 
disturbing force and the inertial reaction increase relative importance. 
the second place, the unknown resistance and virtual-mass effects various 
fittings the pipe line are likely greater importance than the small 
changes indicated Equation (7). 
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printing Paderi’s discussion this paper 
iProceedings Paper 1832 (October, 1958) the fifteen references used that 
were inadvertently omitted. offer them this time for the 


All these branches research offer several developing possibilities, 
attested the authors surveys, which are very valuable either the 
side because reflections the theoretical field. 
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THE HYDRAULIC JUMP ABRUPT 


Closure Walter Moore and Carl Morgan 


WALTER MOORE! A.M. ASCE and CARL A.M. ASCE.— 


The authors wish thank Mr. McPherson and Mr. Forster for their valuable 


discussions. 


Mr. McPherson has compared and integrated the findings with those 


others. asks about the influence the direction movement the tail- 


water level the limiting depths for the various types jumps. brief 
investigation the effect rising and falling tailwater, made the time 


the first experiments, indicated negligible effect this variable provided 


sufficient time was allowed for equilibrium established. result 


the comments Mr. McPherson and the data presented additional ex- 


periments this effect were made the Hydraulic Laboratory the Uni- 
versity Texas. These experiments indicated that sufficient time al- 


lowed the effect the direction tailwater movement was the same 
magnitude the experimental error determining the tailwater depth for 


Froude numbers less than fact for 


for 


and was possible set the tailwater that the 


jump would slowly alternate between the type and the wave type. order 
change from the type the wave type vice-versa moderate time 
was necessary set higher value with rising tailwater than with falling 
tailwater, result similar those observed others. Those results, how- 
ever, were for different range and relative drop height than 


used the authors’ experiments. The equilibrium condition which the 


rising tailwater the change from the type the wave required 
12.3 while for falling tailwater the change occurred 11.4. 
felt that this may have been due the fact that this Froude number the 
top the wave rose higher than the walls the flume. the other hand 
may that this high Froude number the flow pattern such that the 
“hysteresis effect” mentioned McPherson and others present. 

the additional experiments mentioned above, tailwater depths were as- 
certained means piezometers rather than point gage the 


original experiments. The experiments with were made scale 


Proc. Paper 1449, December, 1957, Walter Moore and Carl 
Morgan. 


Chmn., Dept. Civ. Eng., Univ. Texas, Austin, Tex. 
Asst. Prof., Dept. Civ. Eng., Univ Texas, Austin, Tex. 
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twice that used the original experiments. These new tests indicated that 
the values shown Fig. for the maximum jump are slightly high. 
here presented better define the several values required insure 

Mr. McPherson calls attention the disagreement with data presented 
Hsu(4,5) which systematic transition between the types jump (A,B, 
Wave) was found function the Froude No. Actually the investiga- 
tion reported the authors was started explore some longitudinal charac- 
teristics the jump drop assuming the existence the systematic 
transition reported Hsu. was soon found that type jumps could ob- 
tained low values where Hsu’s curves indicated only jump would 
form. Further investigation showed that the full range jump types was 
possible any value (at least above 2.0) controlled only the tail- 
water elevation. The systematic transition function Froude No. was 
not observed the authors. 

From considerations similar those leading figure apparent 
that tailwater control which produced relatively flat curve tailwater 
elevation discharge opposed the steep one shown Fig. would 
result systematic transition from Jump through the wave Jump 
the discharge was decreased. With the depth held constant reduction 
discharge would correspond reduction Fy. This could explain 
systematic change jump type with decreasing reported Hsu, 


However, not clear how the slope the curve tailwater depth 


Froude No. could positive for the type jumps and negative for 
the wave type jump indicated his Fig. Reference 

Mc. McPherson’s Fig. D-2 calls attention again the fact that the band 
shown for Jump Figs. and does not represent the maximum 
values the tailwater. the Results section the paper stated that 
“The conditions under which Jump will form are shown very narrow 
band, but actually Jump will form values greater than the 
upper bound------ 

response Mr. McPherson’s suggestion values are listed 
below for the various forms the jump for which relative bottom velocities 
are presented Figs. 10, 11. 


Mr. Forster effectively emphasizes the importance the type jump 
(type wave, type determining the severity attack the channel 
bottom just downstream from the drop. His examples well illustrate the ef- 
fectiveness Jump and the wave protecting the bottom from high veloci- 
currents. However, the type can result some fairly high bottom 

Figures 10, and the various regions and curves are identified letter 
and number listed the table. can seen that with the wave type and 


3.0 and for 4.0. may seem odd that Froude number 


7 
resulted higher relative bottom velocities than Froude numbers 
and This does not seem unreasonable, however, when the complex flow 


pattern below the drop considered. With large and small Froude 
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FROUDE NO.= 


FIG. FORMS THE HYDRAULIC JUMP 


FUNCTION FROUDE NO. AND RELATIVE 
DOWNSTREAM DEPTH 
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Oo 
3.1 5.4 8.1 
4.3 6.8 9.4 
Oo 
Fig. 
Zo 
Y; =4 Fig. 11 
4.2 5.9 8.80 10.5 
6.5 8.9 11.4 


number the momentum the jet simply not sufficient penetrate down 
the bottom. intermediate Froude number the jet falling 
wave still has sufficient momentum penetrate down the bottom. 
higher Froude numbers the rise tailwater depth sufficient lift the 
off the bottom, thus preventing high velocities from developing the bottom. 
Mr. Forster points out that overdimensioned drop may unconserva- 
tive that Jump may produced with associated high bottom velocities. 
This certainly true for the conditions assumed; that is, that 
constant. This equivalent assuming that the size the drop con- 
trolled the elevation the upstream channel while the downstream channel 
elevation and tailwater depth remain fixed. If, the other hand, the size 
the drop were controlled lowering the bottom elevation the downstream 
channel for suitable distance while the tailwater elevation remains constant 
controlled conditions farther downstream, oversize drop will then 
conservative. Using the same initial conditions Mr. Forster’s example 
wave type jump (bordering Jump B). the drop increased AZ, 
lowering suitable length the downstream channel bottom, the values 
(bordering Jump A). this manner the low bottom velocities would 
served even with oversize drop. 
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DISCUSSION 


Although Figs. 10, and give some very valuable data indicating the 


range velocities expected near the channel bottom, they are not in- 


tended used for precise estimates velocities particular instances. 
Due the complex nature the flow pattern there necessarily considera- 
ble scatter the data. addition the scour potential the flow certainly 


dependent upon the turbulence characteristics the flow well the mean 


velocity. Thus the velocities should thought merely index the 


scour potential for comparisons between the various types jumps. 
Mr. Forster points out his example that Jump produces bottom 
velocity the same order magnitude drowned jump flat apron. 


The peak values the bottom velocities Figs. 10, and are the 


ita 


order which the same that given Mr. Henry (Trans. ASCE, 
1950, 690, Fig. 33) for moderately drowned jump. These relative high 
velocities are not necessarily harmful paved apron provided down- 
stream from the drop for sufficient distance protect the channel bottom. 
For situation which rising discharge and consequent rise tailwater 
would cause increasing submergence the jump its size increases, the 
introduction abrupt drop would beneficial preventing severe drown- 
ing. Instead drowning the jump would shift toward the wave and Jump 
with consequent reduction attack the channel bottom. 
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FLOOD FREQUENCIES DERIVED FROM RAINFALL 


(1958) Journal the Hydraulics Division, Proceedings Paper 1856, there 
appears discussion Clark Proceedings Paper 1451 that the 
authors were unaware when they prepared their closing remarks the 
discussion Messrs. Riggs and Benson (Proceedings Paper 1690, June, 1958). 
The following remarks, therefore, are the closing remarks Mr. Clark’s 
discussion and form addendum the closure published November: 
Generally speaking, the differences between Clark’s opinions and those 
the authors are philosophical well statistical nature. The writer 
criticizes the authors’ use the term “unusual” describing the 1933 
Codorus Creek flood. Clark claims that the most intense 6-hour rainfall 
over square miles the storm producing this flood, i.e., 6.8 inches, 
not sufficient approximation the estimated probable maximum value 
inches called unusual. that measure unusualness? 
Methuselah reputed have lived 900 years. Hence, following Clark’s 
reasoning, would not unusual for someone live 200 years! 
The results recent extensive rainfall-frequency study(1) involving all 
available precipitation records the Middle Atlantic region indicate that the 
shortest average recurrence interval that could reasonably assigned the 
6.0 inch, 10-square mile, 6-hour rainfall the October 1942 storm for either 
Codorus Passage Creek basins would about 160 years. Now, there 
reason why rainfall that intensity should always occur when the sur- 


face layers the soil are, practically speaking, saturated, they were 
the case the October 1942 Passage Creek flood. The same rainfall inten- 


sity could occur under much less favorable flood-producing soil 
Consequently, the shortest average recurrence interval the October 1942 
flood can reasonably expected longer than 160 years, indicating that 
the writer’s flood-frequency curve for Passage Creek probably yields ex- 
cessive peak discharges for the longer recurrence intervals. Examination 
the streamflow records reveals that the curve very likely yields excessive 
discharges for the shorter recurrence intervals also. For example, the curve 
shows two-year peak discharge about 6000 cfs although that value has 
been equaled exceeded only three times the 26-year period, 1932-57 (see 
Basin Table below). 


Proc. Paper 1451, December, 1957, Paulhus and Miller. 
Staff Hydrologist, Weather Bureau, Washington, 
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The shortest average recurrence interval that might reasonably as- 
signed the 7.4 inch, 10-square mile, 6-hour rainfall the 1899 storm men- 
tioned the writer about 800 years. The peak discharge 36,000 cfs 


that Clark estimates could happen Codorus Creek from such rainfall 


under optimum flood-producing conditions would very likely have longer 
average recurrence interval than 800 years. should realized, however, 
that even the average return period such flood were 1000 years, per- 
cent the intervals between occurrences might less than years, and 
percent might less than years.(2) Obviously, the occurrence 

Codorus Creek peak discharge this magnitude within the next few 

years should not astonishing. Nevertheless, the event may still con- 
sidered unusual. 

Now, how about the unusualness the August 1933 maximum peak dis- 
charge 11,200 cfs observed Codorus Creek? This was still the maxi- 
mum observed for the 25-year period ending with water year 1957. According 
the laws probability there 1-percent chance that the average return 
period this flood peak could over 2000 years, and 5-percent chance 
that could about 500 years.(2) Consequently, the 11,200 cfs peak dis- 
charge could unusual event. 

The chief trouble with Clark’s line reasoning that there objec- 
tive way for determining when where stop. For example, the Smethport, 
Penn., storm July 17-18, 1942 (OR 9-23) produced the maximum 10-square 
mile, 6-hour rainfall ever observed the United States, namely, 24.7 inches. 
For some reason, the writer does not mention this storm. Perhaps was 
believed that rainfall that intensity could not occur over either Codorus 
Creek Passage Creek basins. Why not? The probable maximum rainfall 
for both basins estimated somewhat higher than for the storm site. 
Presumably then, rainfall equal greater intensity than for Smethport 
storm could occur could have occurred over these two basins. The fact 
remains, however, that there evidence the occurrence such 
rainfall either these basins within the period our recorded precipita- 
tion streamflow history, and the authors know way base frequency 
analysis what “could happen.” 

Referring Passage Creek, Clark states, “This stream had its floods; 
Codorus Creek did not. least not the period 1933-54. the 26-year 
period 1932-57, Passage Creek experienced three peak discharges exceeding 
the maximum for Codorus Creek. not possible that Codorus Creek 
less prone flooding than Passage Creek? the other hand, not pos- 
sible that Passage Creek, through chance distribution, has received more 
than its share floods for this 26-year period? Table suggests that the 
latter definitely possible. Passage Creek the only basin for which the 
three highest annual flood peaks are times the fourth highest. fre- 
quency analysis (Gumbel method adjusted for partial duration series) its 
annual floods for the 26-year period, 1932-57 yields the following: 


Peak discharge (1000 cfs) 4.5 8.5 11.5 15.5 18.6 21.5 


These data differ from Clark’s frequency curve, presumably based re- 
gional considerations, less than percent except for return periods 
shorter than about years. Comparison the statistical arrays annual 
floods for the seven basins Table indicates that the Passage Creek record 
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Table 


Rank 


COW ON DU FWP 


Chester Creek near Chester, Penn., sq. mi., 1932-57 
Codorus Creek Spring Grove, Penn., sq. mi., 1933-57 
Linganore Creek near Frederick, mi., 1935-57 


Seneca Creek Dawsonville, Md., sq. mi., 1931-57 


DISCUSSION 


STATISTICAL ARRAY ANNUAL PEAK DISCHARGES (100 CFS) 


FOR SEVEN BASINS SQUARE MILES 


Passage Creek Buckton, Va., sq. mi., 1932-57 


Creek near Pipersville, Penn., sq. mi., 1936-57 
Tye River near Lovingston, Va., sq. mi., 1939-57 


provides about the most unrepresentative sampling for the group. Can Linsl 
Clark’s curve then accepted representative for this group basins? 
The authors were unable determine the plotting formula used the “Seas 
criteria for construction the curve the writer’s frequency plot. Neither 105th 
could determined why the frequency curve ascribed the authors was 
plotted from their Table For the Codorus Creek flood peaks plotted, the 1956. 


curve should have been drawn from data the authors’ Table Neither 
set data recommended representative for Codorus Creek. For the 
situation under discussion, the authors recommend the use all 
streamflow records for the basin supplemented much possible syn- 
thesized flood peaks. However, the synthesized flood peaks for the full 


period rainfall record should yield good results. Table compares these 
results with those from the observed record and shows that the difference 
appreciable. The data were taken from Tables and the authors’ 
original paper. 

Table FLOOD FREQUENCIES (100 CFS) CODORUS CREEK 

(By Gumbel Analysis Not Adjusted for Partial Duration Series) 


Return Period (years) 


view Clark’s questioning the reliability flood frequencies derived 
from something other than floods and from data outside the watershed, the 
authors were astonished his reference rainfall data outside the subject 
basins indication the possible frequency extreme values ob- 
served peak discharges. Furthermore, the manner which Clark uses 
streamflow and/or rainfall data within large region establish flood fre- 
quencies for small basin suggests tendency confuse the frequency 
occurrence extreme event within the region with the probability 
its occurrence the small basin. 

conclusion, next streamflow measurements, rainfall data are 
the best index flood frequencies. The idea regional procedures sound, 
but there reason why such procedures based the relatively short 
records streamflow generally available should reliable for the longer 
recurrence intervals. The authors’ procedure could used increase 
synthetically the length streamflow records for individual basins pro- 
vide more reliable base for regional procedures that would yield much more 


dependable results, especially for the longer recurrence intervals, than 
now possible. 


“Rainfall Intensity-Frequency Regime, Part The Middle Atlantic 
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Region,” Weather Bureau Tech. Paper 29, press. 
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TURBULENCE CHARACTERISTICS THE HYDRAULIC JUMP* 


Closure Hunter Rouse, Siao, 
and Nagaratnam 


ments upon this paper were particular interest that they reflected the 
views investigators five the principal hydraulics laboratories now en- 
gaged research turbulent flow. Significant points already considered 
during the preparation the paper were emphasized, and others not con- 
sidered least not given their proper attention were very pertinently 
noted. always the case, the net result marked improvement over 
the original contribution. 

The writers recommend that everyone who has access glass-walled 
flume carry out the experiment described Mr. Silberman artifically 
varying the velocity distribution, necessary, the application coarse 
bed roughness. The formation standing eddy level floor con- 
vincing demonstration basic boundary-layer principle one could 
hope find. Other hydraulic distortions the jump proper have been dis- 
cussed detail elsewhere, but they might well mentioned here this 
connection: those due variation channel cross section, particularly 
through the trapezoidal range; variation bed slope, chute, bed 
apron; and variation wall alignment curvature, flared warped 
chute spillways. Each these will have influence upon the jump charac- 


teristics just important its way that bed roughness but just in- 


essential the basic jump phenomenon. conducting the investigation, the 
writers attempted eliminate every possible extraneous factor including 
initial velocity distribution that full attention could given the 
mechanism turbulent expansion. Surely the introduction sonic throat 
and shock wave, suggested Mr. Silverman, would not have been ef- 
fective step this direction. 

Messrs. Peterka and Bradley introduced welcome note into the literature 
their immediate recognition the usefulness academic findings 
clarifying the practical aspects related problem. The question jump 
length has long been matter controversy, part for lack agreement 
nomenclature and part because the jump actually never ends but only 


Proc. Paper 1528, February, 1957, Hunter Rouse, Siao and 
Nagaratnam. 


Hydr. Engr., Inst. Hydr. Research, Academia Sinica, Peking, China. 
Engr., Harza Engr. Co., Chicago, 


Director, Inst. Hydr. Research, State Univ. Iowa, Iowa City, Iowa. 
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approaches asymptotically state depth, velocity, and turbulence equilibri- 
um. Designation end thus purely arbitrary. Messrs. Peterka and 
Bradley have realistically adopted criterion based upon bed velocity; 
Bakhmeteff and Matzke logically based theirs depth; the writers con- 
veniently terminated their boundary curves where the effect further curva- 
ture seemed inappreciable. Were residual turbulence (the academic 
choice) determine the nominal end, its position would obviously depend 
upon the percentage arbitrarily chosen for the permissible turbulence re- 
mainder. Since this remainder continues appreciable comparison 
with that the subsequent uniform flow for some distance downstream, such 
definition would produce far the greatest nominal length. The discussors 
are right their belief that length roller and length rapid expansion 
would approach equality the Froude number became very large; the 
residual turbulence were related the overall loss, all definitions would 
then agree the (unfortunately, purely academic) limit. 

Mr. Robertson, select minor comment from his discussion this 
point, speaks normal “proper” hydraulic jump one which occurs 
established open-channel flow. This recalls, course, the discussion Mr. 
Silberman. Messrs. Peterka and Bradley, the other hand, may well think 
proper jump one controlled the baffle piers foreshortened 
apron. “Established” flow actually more definite than the design 
pier and apron, for can vary through wide limits with channel roughness 
and shape. For this reason Mr. Harleman’s detailed momentum analysis 
the jump, and his application the resulting equation the three conditions 
described the paper, represent supplementary material basic import. 
The writers fear, however, that readers will ignore the equation and assume 
the plotted curve indicate general solution. The fact should hence em- 
phasized that even established flow there will different deviation from 
the elementary equation for every Reynolds number, relative roughness, and 
channel shape. 

return Mr. Robertson’s discussion, the writers can only report that 
there was indication bottom separation either flume duct for the 
condition minimal boundary-layer growth that was sought experimentally. 
far the upper boundary effect was concerned, the vicinity the non- 
uniformity was surely not one shear; whereas more refined method 
measurement would undoubtedly have revealed positive shear above the 
roller, correctly expected Mr. Robertson, comparison its probably 
magnitude with that the region high velocity gradient below would not 
lead one anticipate much change from the conditions shown. The effect 
the boundary upon large-scale fluctuations, the other hand, was admittedly 
considerable, and some indication its magnitude undoubtedly found 
Mr. Hubbard’s measurements. 

Free-turbulence literature now distinguishes among three different re- 
gions the turbulence spectrum: the fine-scale region, containing the eddies 
which most the energy the turbulence concentrated; and the large- 
scale region intermittency described Mr. Robertson. this sense the 
large eddies which are obvious part the true hydraulic jump represent 
neither appreciable concentration energy nor important source 
dissipation. They can, however, play considerable role the process 
mixing. the juncture between the roller and the main stream the intermit- 
tent formation these large eddies (often distinguishable jump photo- 
graphs) much like that observed the edges jets, wakes, and boundary 


no 


layers. The free surface, which the more noticeable zone intermittent 
fluctuation and that which Mr. Robertson refers, has some the same 
characteristics. However, noted the discussion, also has two inter- 
related features which distinguish it: abrupt discontinuity the liquid 
medium itself, and the stabilizing restoring action gravity. The air 
model surely duplicated the essential freedom motion the lower zone 
intermittency; what extent the fluctuations this zone are actually con- 
trolled the free-surface fluctuations remains determined. 

Mr. Robertson correctly implied that the plotted dissipation functions were 
partially assumed rather than completely measured nature, view 
the representation the nine velocity derivatives single one through the 
relationships isotropy. Whether this assumption requires only high Rey- 
nolds number closely valid, the components even then would deviate 
different degrees different directions and different zones flow, must 
also remain future investigations for conclusive answer. Surely the 
present measurements not indicate the existence much variation 
validity from section section, else the production-convection -dissipation 
sum would have departed more radically than did from zero value. 
should also remarked this regard that such indications small-scale 
isotropy exist date have been found under conditions free turbulence, 
whereas deviations therefrom have all been noted the vicinity solid 
boundaries. 

Measurements the type proposed implicitly Mr. Robertson and ex- 
plicitly Mr. Hubbard must eventually made, although the writers shud- 
der think the painstaking labor involved determining many succes- 
sive sections under conditions which even the depth itself very hard 
define what date has been measured only for representative sections 
flows which are either uniform similar from section section. Mr. Hub- 
bard apparently has hand equipment that will accomplish some part the 
measurements required least safe distance below the freely 
fluctuating surface and has demonstrated knack filling other needs 
they arise. Still found are group expert laboratory technicians 
with considerable time their disposal, and unusually patient and 
sympathetic sponsor! 
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SNOWMELT 


KENNETH A.M. ASCE.—Snowmelt forecasting usually 
considered being limited the mountainous regions the country. 
this paper the authors have shown that can also utilized non- 
mountainous areas. 

The method presented the authors appears sound approach 
the problem snowmelt forecasting, and does present procedure for 
developing similar relationships for other non-mountainous areas. 

The writer has checked the reliability the relationship 


Where 


applying three years records much the same manner river- 
level forecaster would use it. The three years chosen were 1953, 1954, and 
1955, which were subsequent the authors’ period study. 

The quarter-day dry bulb temperature, dew point temperature, wind speed, 
and sky cover for the Wausau, Wisconsin, station were used calculate the 
quarter-day melt potentials. The melt potentials were summation the 
convection, condensation, and radiation melt potentials, was stated the 
authors. 

Precipitation was also determined quarter-day basis throughout each 
the three spring seasons. The Medford and Marshfield, Wisconsin, sta- 
tions were used determine the precipitation. These are the two stations 
used the authors. The rainfall was added directly the total quarter-day 
melt potentials. 

Daily discharges the Big Eau Pleine River given the 
Geological Survey were corrected according the “ice index theory,” given 
the paper, for the period ice effect. The beginning the melt period 
and the total inches runoff were determined from these daily discharges. 

The meteorlogical data was then used river-level forecaster would 
use the data determine future flows. The forecasted flows were then 
checked against the actual daily river discharges. Due the unusually small 
snow accumulation 1954, however, was necessary resort the 
authors’ Table forecast the snowmelt hydrograph for 1954. 


Proc. Paper 1834, November, 1958, Zoller and Lenz. 
Cons. Engr., Boulder, Colo. 
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The forecasted hydrographs are compared the recorded hydrographs 
Figures and for 1953, 1954, and 1955. 1955 the forecasted and re- 
corded hydrographs are quite close for the first five days the snowmelt, 
and then they diverge. The divergence due the fact that the snow the 
ground had completely melted. March and 22, though, there was new 
snowfall which required ripening, and which then ran off March and 
indicated the hydrograph. The forecasted hydrographs were based 
the runoff reaching the gaging station the day after the melt potential oc- 
curred. This accordance with the authors’ statement that there 
13-hour lag runoff for the basin and also that most the melt potential 
occurs the second period each day. 

The forecasted and recorded runoff for the three years was put into mass 
curves and also compared. These are shown Figures The hori- 
zontal portions the forecasted mass curves are due there being melt 
potential for those specific days. This illustrated for 1955 from March 
through March 19. 

Figure the authors’ paper the relationship developed presented 
graphically and compared with several years reliable data. Figures and 
this discussion give the relationship determined for 1953 and 1955. 

shown that the equation which was developed fits well with the recorded 
runoff for these two years. 1954 was not compared because the very small 
amount snow the ground caused distortion the relationship. 

The summary the melt conditions for the three years studied for this 
discussion follows: 


Length Rain During Inches Runoff 
Year Melt Period Melt Period Melt Period During Melt Period 
1953 March days 0.12 inches 4.17 inches 
March 

1954 March days 0.58 inches 0.76 inches 
March 

1955 March days 1.15 inches 3.77 inches 
April 


The Climatological Data bulletins report that 1953 there were 
inches snow the ground before melting began. 1954 there was only 
about one inch snow, resulting only very small amount runoff. The 
1955 bulletin lists inches snow before melting, and then March 

From the writer’s investigation the three years subsequent the 
authors’ period study, the following conclusions were drawn regard 
the Zoller-Lenz method forecasting snowmelt runoff. 


The forecast the daily discharges for 1953, 1954, and 1955 correlated 
satisfactorily with the recorded discharges the Big Eau Pleine River. 


The relationship developed does not hold true for the runoff from un- 
usually small snow accumulations such 1954. 


With small drainage basin the size the Big Eau Pleine Basin, the 
runoff forecasts can only made for one day advance. The Weather 
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Bureau’s five-day forecasts cannot alleviate this any great extent because 
the amount detailed meteorlogical data needed for the authors’ method. 


The determination when the snow will attain ripeness and thus when 
additional melt potential will cause runoff problem which has not been 
solved this paper. small drainage basins such the Big Eau Pleine 
Basin, the forecaster would probably miss the first day’s snowmelt runoff. 
would, however, alerted the ripeness this first day’s runoff and 
then begin his forecasts for the subsequent days. 


The “ice index theory” appears sound approach the problem 
determining discharges during the rising stage the hydrograph. 


summary, this paper snowmelt runoff valuable contribution 
the published information hydrology. The fact that this study was based 
non-mountainous region indicates that the interest snowmelt runoff 
spreading beyond the areas where comprises nearly the entire supply 
water. 
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DIVISION ACTIVITIES 
HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 


February, 1959 


PURPOSE THE HYDRAULICS DIVISION 
(Quoted from the Official Register) 


“The advancement and dissemination knowledge relating the occur- 
rence water nature and its behavior structures, water courses, and 
underground. 

“In particular the field the Hydraulics Division shall embrace meteorolo- 
and hydrology they affect the engineer, fluid mechanics engineering 
usage, and applied hydraulics branch engineering science which fur- 
nishes the basis for hydraulic design and for the practical use water the 
different specialized branches hydraulic engineering.” 


MEMBERS, ASCE HYDRAULICS DIVISION EXECUTIVE COMMITTEE, 
1958-1959 


Term Expires 


Professor Carl Kindsvater, Chairman 1960 
George Institute Technology 
Atlanta, Georgia 


Dr. Arthur Ippen, Vice-Chairman 1961 
Professor Hydraulics 

Massachusetts Institute Technology 

Cambridge, Massachusetts 


Mr. Harold Martin 1959 
Bureau Reclamation 

Building 

Denver Federal Center 
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Note: No. 1959-9 part the copyrighted Journal the Hydraulics Division, 
Proceedings the American Society Civil Engineers, Vol. 85, February, 
1959. 


Copyright 1959 the American Society Civil Engineers. 
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Mr. Earl O’Brien, Contact Member from ASCE 1958-59 
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O’Brien and Gere 

400 East Genesee Street 
Syracuse New York 


COMMITTEE CHAIRMEN, ASCE HYDRAULICS DIVISION, 1958-59 


Committee Hydrology: 


Committee Hydraulic Structures: 


Committee Flood Control: 


Committee Sedimentation: 


Committee Hydromechanics: 


Committee Tidal Hydraulics: 


Committee Publications: 


Committee Research: 


Mr. Charles McDonald 
Department the Interior 
Geological Survey 
Washington 25, 


Mr. Glenn Hands 
Box 7088 
Kansas City, Missouri 


Mr. Joseph Perrey 
1330 Michigan Street 
Indianapolis, Indiana 


Mr. Albert Gildea 
Los Angeles, California 


Professor John Ripken 

University Minnesota 

St. Anthony Falls Hydraulic 
Laboratory 

Minneapolis 14, Minnesota 


Mr. Lindsay Disney 
Coast and Geodetic Survey 
Washington 25, 


Mr. James Ball 
Bureau Reclamation 
Building 

Denver Federal Center 
Denver Colorado 


Professor James Robertson 
125 Talbot Laboratory 
University Illinois 
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Committee Standards: Mr. Rolland Carter 
Geological Survey 
Department the Interior 
Washington 25, 


Committee Stevens Award: Professor David Appel 


University Kansas 
Lawrence, Kansas 


Committee Hilgard Prize: Mr. Albert Fry 


Tennessee Valley Authority 
Knoxville, Tennessee 


MEETING THE EXECUTIVE COMMITTEE, HYDRAULICS DIVISION 


two-day meeting the Executive Committee was held St. Louis 
November 3-4. Those attending were: 


November and 


Professor Carl Kindsvater, Chairman 
Dr. Arthur Ippen, Vice-Chairman 
Mr. Harold Martin 

Mr. Maurice Dickinson 

Mr. Joseph Tiffany, Secretary 


November 


Mr. Glenn Hands, Chairman, 
Committee Hydraulic Structures 


The following items have been extracted from the minutes the meeting: 


HYDRAULICS DIVISION BUDGET. Approval has been receiver from 
ASCE Headquarters for the following 1958-59 budget: 


Technical Committee Activities $11,200 
Administrative Acitivities 2,300 
TOTAL $13,500 


NEWSLETTER. Consideration was given headquarters request that 
the Newsletter schedule conform with more frequent mailing schedule the 


Journal; however, the concensus was that the Newsletter should still pub- 
lished once every two months. 


COMMITTEE THE KARL EMIL HILGARD HYDRAULIC PRIZE. 
change the wording the purpose this committee, submitted Mr. 


Tiffany, was approved. The revised statement will appear the 1959 Official 
Register. 


THEME FOR BOSTON MEETING. The Committee concurred with the sug- 
gestion Mr. Vincent, Structural Division, ASCE, regarding the adoption 
the theme, “Basic Research and Structural Engineering,” for our portion 
the technical program the Boston meeting. 


FUNDS FOR BASIC RESEARCH. The Committee unanimously approved the 
resolution the Structural Division, ASCE, proposed for presentation the 
Technical Procedure Conference 1959, support the establishment 
ASCE funds for basic research engineering sciences. 
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NEXT MEETINGS THE EXECUTIVE COMMITTEE. Two additional 
meetings the Executive Committee are planned for the remainder the 
Society year. One meeting will held early 1959 plan the Hydraulics 
Division organizational changes required conform with recommendations 
the Water Resources Task Committee (see below, also October 1958 
Newsletter). The other meeting planned for July the time the Ft. 


Collins conference (see below), and will result next year’s committee 
appointments. 


ASCE CLEVELAND CONVENTION 
May 4-8, 1959 


TENTATIVE HYDRAULICS DIVISION SESSIONS 


(1) Sponsored the Flood Control Committee 


Presiding: Harold Martin, Member, Executive Committee, and 


Joseph Perrey, Chairman, Flood Control Committee, 
Hydraulics Division 


CONSERVANCY DISTRICTS FLOOD CONTROL AGENCIES, 
Chambers, consulting engineer, Dayton, Ohio. 


OPERATION THE MIAMI CONSERVANCY DISTRICT FOR FLOOD 


CONTROL, Max Mitchell, chief engineer, Miami Conservancy District, 


Dayton, Ohio. 


FLOOD CONTROL AND OTHER BENEFITS MUSKINGUM CON- 
SERVANCY DISTRICT, Tom Shuler, chief engineer, Muskingum 
Watershed Conservancy District, New Philadelphia, Ohio. 


(2) Sponsored the Hydraulic Structures Committee 


Presiding: Harold Martin, Member, Executive Committee, and 
Glenn Hands, Chairman, Hydraulic Structures 
Committee, Hydraulics Division. 


USE SUBMERGIBLE GATES NAVIGATION STRUCTURES, 
Abbott, and Moors, Ohio River Division, Corps 
Army, Cincinnati, Ohio. 


MODEL TESTS SUBMERGIBLE TAINTER GATES, Thomas Murphy, 
Chief, Hydraulic Structures Section, Waterways Experiment Station, 


Corps Engineers, Army, Vicksburg, Mississippi. 


RELATIVE PERFORMANCE CURRENT METERS GAGING THE 
DISCHARGE THE OUTFALL RIVERS THE GREAT LAKES, 
Wayne Townsend, and Frank Blust, Hydraulics and Hydrology Branch 
Lake Survey, Corps Engineers, Army, Detroit, Michigan. 


(3) Sponsored the Hydromechanics Committee 


Presiding: Harold Martin, Member, Executive Committee, and 


John Ripken, Chairman, Hydromechanics Committee, 
Hydraulics Division 
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FLOW CONDITIONS THE ENTRANCE PIPE, Harold Henry, 
Assistant Professor Civil Engineering, Michigan State University, 
East Lansing, Michigan. 


BEHAVIOR BUBBLES VIBRATING LIQUID, Richard Skalak, 
Civil Engineering, Columbia University, New York, New York. 


LABORATORY DEVELOPMENT THE CONTRA COSTA CULVERT 
OUTFALL ENERGY DISSIPATOR, Russell Keim, Instructor, Depart- 


ment Theoretical and Applied Mechanics, University Illinois, Urbana, 
Illinois. 


Sponsored the Sedimentation Committee 


Presiding: Harold Martin, Member, Executive Committee, and 
Fred Larson, Past-Chairman, Sedimentation Committee, 
Hydraulics Division. 


THE FUNCTIONS AND OBJECTIVES THE SEDIMENTATION SUB- 
COMMITTEE THE INTER-AGENCY COMMITTEE WATER 
RESOURCES, Warren Murphy, Director, Division Flood Prevention 
and River Basin Programs, Forest Service, Washington, 


OBSERVATIONS FLOOD FLOW EFFECTS CHANNEL BOUN- 
DARIES, Donald Parsons, Hydraulic Engineer, Agricultural Research 
Service, Department Agriculture, East Aurora, New York. 


THE NEW ARS SEDIMENTATION RESEARCH LABORATORY, OXFORD, 
MISSISSIPPI—PROBLEMS AND OBJECTIVES, Russell Woodburn, 
Project Supervisor, Agricultural Research Service, Department 
Agriculture, University, Mississippi. 


HYDRAULICS DIVISION CONFERENCE FORT COLLINS 


The eighth annual conference the Hydraulics Division will held 
Fort Collins, Colorado, July 1-3. Co-hosts are the Colorado Section and 
Colorado State University. Six half-day technical sessions will include 
papers various hydraulic subjects, and tour through the Hydraulics 
Laboratory Colorado State University. Preceding the conference visit 
the Bureau Reclamation Hydraulics Laboratory Denver planned for 
June followed guided tour the Colorado-Big Thompson June 30. 
The Colorado-Big Thompson project, built the Bureau Reclamation, 
one the most extensive irrigation and power projects Colorado. 

Housing arrangements for individuals and families new, modern Uni- 
versity dormitories will available reasonable rates. Social events are 
planned and interesting schedule being arranged for the ladies and 
children. Numerous recreational attractions, including Rocky Mountain 
National Park, are within hour’s drive Fort Collins. Plan come and 
enjoy Colorado vacation trip. you wish reserve mountain cabin for 
family vacation, write Housing Committee, Hydraulics Conference, Civil 
Engineering Department, Colorado State University, Fort Collins, Colorado. 
Reservations for cabins should made early possible because sum- 
mer tourist demand. 


More detailed information planned for the April issue the Newsletter. 
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OCTOBER SESSIONS ASCE BOARD DIRECTION 


During the October sessions the ASCE Board Direction number 
actions were taken particular interest members the Hydraulics 
Division, follows: 


TASK COMMITTEE WATER RESOURCES, ASCE. The Board approved 
the recommendations made the report the Task Committee and acted 
set the coordinating Committee Water Resources recommended 
the report (see October 1958 Newsletter). 


COMMITTEE TECHNICAL PROCEDURE, ASCE. Appointments the 
Committee for 1959 confirmed the Board Direction included Profes- 
sor Carl Kindsvater, Chairman, the Hydraulics Division Executive Com- 
mittee. The Board also acted approve meeting 1959 the Committee 
Technical Procedure, Chicago other central location, with attendance 
the coordinating committee Water Resources. 


TECHNICAL DIVISION CONTACT MEMBERS. Appointments contact 


members for 1958-59 were confirmed the Board Direction October. 


Contact member for the Hydraulics Division Mr. Earl Direc- 
tor representing District 


COMMITTEE SOCIETY FELLOWSHIPS, SCHOLARSHIPS, GRANTS 
AND REQUESTS, ASCE. This committee was created action the Board 
Direction, 13-14 October, amendment the Bylaws (Art. VII, Sect. 
3(f). Members this new committee are: 


Professor John Vennard, Chairman Oct. 1959 
Dr. Lorenz Straub Oct. 1960 
Dr. Martin Mason Oct. 1961 
Dr. George Hickox Oct. 1962 
Mr. Samuel Morris Oct. 1963 


matter interest, four these men are members our division, 
and two are past chairmen the Executive Committee. 


FELLOWSHIPS AND GRANTS 


Applications for the FREEMAN FELLOWSHIP are now being accepted for 
the academic year 1959-1960, and can submitted any time before March 
1959. Younger members ASCE ASME are eligible for this award. 

new fellowship grant the field CIVIL ENGINEERING RESEARCH has 
just been announced ASCE. The grant the amount $5,000 made an- 
nually and for the purpose aiding the creation new knowledge for 
the benefit and advancement the science and profession civil engineering. 
Applicants must members ASCE, citizens, and have been gradu- 
ated from accredited curriculum. 

Correspondence concerning the above fellowships should addressed 
the Executive Secretary, ASCE, West 39th Street, New York 18, New York. 


The Eighth Congress the International Association for Hydraulic Re- 
search will held the Windsor Hotel Montreal, Canada, during the week 
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beginning August 24, 1959. 

All members and those interested the work the Association are in- 
vited attend and contribute the meeting, both presentation papers 
and participation the discussion. 

The following four subjects have been selected for study the technical 
sessions the General Assembly the Lisbon Convention 1957. 


Hydraulics Gates and Valves (including vibration, forces, cavitation, 
ice). 


Fundamental Hydraulics Ship Locks. 


Density Currents (including effects reservoir sedimentation, hydro- 
electric intakes, fishing industry, etc.) 


Air Entrainment and Air Vents. 


Two seminars for more informal study will conducted concurrently with 
the technical sessions the following subjects: 


Ice Problems Hydraulic Structures. 


Transportation Material Water. (It intended that this seminar 
deal primarily with sediments but transport other material may 
included.) 


addition the technical sessions, excursions are being planned the 
St. Lawrence Seaway structures, some the large hydro-electric projects 
the region and hydraulic laboratories that have contributed the de- 
velopment both. this connection noted that the St. Lawrence 


Seaway Project will officially opened 1959 and will thereafter provide 
passage for deep draft ocean-going vessels the heart the North American 
Continent. 

program being planned for the ladies that will enable them become 
better acquainted with Montreal, the second largest French speaking city 
the world. 


The official languages the Congress will French and English. 


TECHNICAL PAPERS. Those intending submit papers any the 
four official subjects should advise the Canadian Committee January 
15, 1959, giving the title and brief description the contents the paper 
each case. 

The full texts the papers must the hands the Canadian Com- 
mittee not later than March 1959, order guarantee publication and 
distribution advance the Congress. 

Each paper include summaries not exceeding 300 words both 
French and English. 

Papers that have been translated into one the official languages should 
accompanied copy the original language. 

Every author attending the Congress will invited present outline 
his paper resume not exceeding 500 words length. 

Notice intention participate either the seminars should re- 
ceived the Canadian Committee not later than March 1959. Copies 
prepared contributions either seminar should received not later than 
June 1959, ensure being included the Agenda. 
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PROCEEDINGS. expected the proceedings the Congress will 
published. These proceedings will include: 


Papers the four official subjects. 

Discussion the four official subjects. 

Minutes the Seminar. 

Free papers accepted the Executive Committee the Council. 


FEES. The fee for the Congress $20.00 for participants and $10.00 
(Canadian Currency) for the ladies, exclusive costs excursions. the 
case non-members, the fee will $25.00 for participants, which include 
membership, for those who are eligible, one year the Association. (The 
fees for ladies remaining the same.) The range prices per day for rooms 


Single room with bath 8.00 $13.00 
Double room with bath $11.00 $17.00 


The cost meals extra and may expected vary from $1.00 
$3.00, varying with the tastes the individual. 
Travel and hotel reservations may made through Thos. Cook and Sons. 


COMMUNICATIONS. Additional circulars will issued due course 
giving details the program and instructions concerning the preparation 
papers, registration, and other matters. All communications regarding at- 
tendance, participation, submission papers, etc., should addressed to: 
Mr. Leo Roy, c/o Quebec Hydro-Electric Commission, 107 Graig Street West, 
Montreal P.Q., Canada. 


MONTHLY SCIENTIFIC PUBLICATION INTEREST 
‘TO HYDRAULICS DIVISION MEMBERS 


will interest our Division members know that the American 
Geophysical Union combining the scientific and technological material 
its bimonthly TRANSACTIONS with the quarterly JOURNAL GEOPHYSI- 
CAL RESEARCH and will issue the JOURNAL GEOPHYSICAL RESEARCH 
monthly starting 1959. This will contain much material bearing 
Hydrology interest our members. 

The subscription rate for the new JOURNAL $16.00 for the calendar 
year. Memberships, however, are invited. Membership $10.00 per calen- 
dar year and includes subscription the JOURNAL. Any who are inter- 
ested should write the office the American Geophysical Union, 1515 
Massachusetts Avenue, W., Washington 


FOR YOUR CALENDAR 
ASCE Meetings 
May 4-8, 1959 ASCE, Cleveland Convention 
July 1-3, 1959 Hydraulics Division, Fort Collins Conference 
October 19-23, 1959 ASCE, Washington, C., Convention 
March 7-11, 1960 ASCE, New Orleans Convention 


June 19-23, 1960 ASCE, Reno Convention 
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ASCE Hydraulics Division 1959-9--9 


October 9-13, 1960 ASCE, Boston Convention 
April 10-15, 1961 ASCE, Phoenix Convention 
October 16-20, 1961 ASCE, New York Convention 
February 1962 ASCE, Houston Convention 
May 1962 ASCE, Omaha Convention 
October 15-19, 1962 ASCE, Detroit Convention 


Non-ASCE Meetings 


February 22-28, 1959 NATIONAL ENGINEERS WEEK 

April 13-15, 1959 ASME Hydraulic Conference, University 
Michigan 

June 15-19, 1959 American Society for Engineering Education, 
Pittsburgh, Pennsylvania 

June 17-20, 1959 National Society Professional Engineers, 
New York City 

August 1959 International Association for Hydraulic 


kesearch, Montreal, Canada 
Deadline dates for Newsletter contributions: June issue April 15. 


THE HYDRAULICS DIVISION NEWSLETTER 


Your editor sincerely appreciates the cooperation the various members 
furnishing timely and pertinent information for the Newsletter. You are 
urged continue use the Division Newsletter for announcements, in- 
quiries, personnel news, committee reports, surveys and other items 
interest Division members. 


GUY ARBUTHNOT, JR. 
Box 631 

Vicksburg, Mississippi 
Newsletter Editor 
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